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PREFACE 


The  study  reported  herein  was  authorized  by  the  Office,  Chief  of 
Engineers  (OCE),  U.  S.  Army,  on  l6  July  1975,  at  the  request  of  the 
U.  S.  Army  Engineer  District,  Sacramento  (SPK). 

This  study  pertains  to  the  900-MW  power  project  and  is  Report  1 
of  a series.  At  the  time  this  report  was  prepared,  another  study  was 
being  made  of  the  same  project  with  a greater  power-producing  capa- 
bility (2250  MW),  The  report  of  that  study  will  be  published  separately 
as  Report  2. 

The  Investigation  was  conducted  during  the  period  August  1975  to 
October  1976  in  the  Hydraulics  Laboratory  of  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  the  direction  of  Messrs.  H.  B. 
Simmons,  Chief  of  the  Hydraulics  Laboratory,  and  J.  L.  Grace,  Jr.,  Chief 
of  the  Structures  Division  and  Reservoir  Water  Quality  Branch  (Physical). 
The  study  was  conducted  by  Messrs.  M.  S.  Dortch,  D.  G.  Fontane, 

B.  Loftis,  and  C.  H.  Tate,  Jr.,  with  assistance  from  Mr.  J.  H.  Riley. 

This  report  was  prepared  by  Messrs.  Fontane,  Dortch,  Tate,  and  Loftis 
and  reviewed  by  Mr.  Grace. 

Representatives  of  WES  and  SPK  met  at  SPK  during  November  1975  to 
discuss  the  scope,  objectives,  and  approach  of  the  study.  During 
July  1976,  representatives  of  WES,  SPK,  and  the  South  Pacific  Division 
(SPD)  met  to  redefine  the  scope  and  objectives  of  the  study.  Mr.  Dennis 
Huff  of  SPK  was  detailed  to  WES  during  September  1976  to  become  familiar 
with  the  mathematical  model  and  to  assist  with  the  numerical  simulations. 

Directors  of  WES  during  this  study  were  COL  G.  H.  Hilt,  CE,  and 
COL  J.  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply 

feet 

miles  (U.  S.  statute) 
square  feet 

square  miles  (U.  S.  statute) 

acres 

acre-feet 

cubic  feet  per  secoiid 
British  thermal  unit 
degrees  Fahrenheit 


■ 

O.30U8 

1.609 

0.092903 

2.58999 

h.OhS 
1233. ij82 

0.02832 

1055.056 

5/9 


To  Obtain 

metres 

kilometres 

square  metres 

square  kilometres 

square  kilometres 

cubic  metres 

cubic  metres  per  second 

Joules 

degrees  Celsius  or  Kelvins* 


fvlARYSVILLE  LAKE  HYDR0THERJ4AL  STUDY 


900-r«IW  PROJECT 

Hydraulic  and  Mathematical  Model  Invest ip:at ion 


PART  I : INTRODUCTION 

Purpose  and  Scope  of  Study 

1.  This  study  was  conducted  to  determine  the  ability  of  the 
Marysville  Lake  project  to  satisfy  downstream  water  temperature  ob- 
jectives. The  study  required  (a)  understanding  and  description  of  the 
large-scale  hydrodynamic  phenomena  within  the  project,  (b)  development 
of  temperature  profiles  within  Marysville  Lake,  and  (c)  determination 
of  temperatures  released  from  Marysville  Lake,  from  the  afterbay,  and 
routed  to  the  mouth  of  the  Yuba  River. 

2.  A physical  hydrodynamic  model  and  a numerical  simulation  model 
were  used  in  the  study.  The  physical  model  provided  information  rela- 
tive to  the  hydrodynamic  response  of  the  prototype  for  various  oper- 
ational conditions.  The  numerical  model  provided  the  capability  for 
assessing  for  year-long  periods  the  effect  of  historical  meteorologic 
and  hydrologic  data  and  various  operating  conditions  on  temperature 
regimes  within  and  downstream  of  the  project. 

Project  Description 

3.  The  Marysville  Lake  project,  authorized  by  Congress  in  the 
Flood  Control  Act  of  7 November  1966  (Public  Law  89-T89)  and  modified 

by  the  Water  Resources  Development  Act  of  1976  (Public  Law  9^-587),  will 
provide  flood  control  for  downstream  areas,  hydroelectric  power,  water 
for  irrigation,  recreation,  and  fishery  enhancement.  The  project  will 
be  built  and  operated  by  the  U.  S.  Army  Corps  of  Engineers.  Irrigation 
and  power  functions  of  the  project  will  be  integrated  into  the  Central 
Valley  Project  of  the  U.  G.  Bureau  of  Reclamation. 
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It.  The  lake  will  be  on  the  Yuba  River  about  17  miles*  upstream 
from  the  city  of  Marysville  (Figure  1)  in  north-central  California.  Tiie 
project  will  consist  of  a multipurpose  pool  impounded  by  main  d!ur;s  on  the 
Yuba  River  just  above  the  Parks  Bar  Bridge  and  on  the  adjacent  Dry  Creek 
and  an  afterbay  dam  and  pool  for  regulation  of  downstream  flows.  The 
afterbay  will  also  be  used  to  store  water  to  be  pumped  back  to  the  lake. 

5.  The  project  was  relocated  to  the  Parks  Bar  site  from  the  author- 
ized Browns  Valley  site  following  a review  of  alternatives  in  light  of 
the  recent  energy  crisis.  Studies  indicated  that  the  Parks  Bar  site  had 
greater  potential  for  hydroelectric  power,  including  poraped-storage. 

Also,  the  Parks  Bar  site  was  found  to  be  .more  environmentally  and 
socially  acceptable. 

6.  The  dam  on  the  Yuba  River  will  rise  368  ft  above  ttie  streambed 
to  el  570**  and  consist  of  a concrete,  gravity-center  section  with  earth- 
fill  abutments.  The  dam  on  Dry  Creek  will  rise  328  ft  above  the  stream- 
bed  and  consist  entirely  of  earthfill.  Both  dams  will  be  7000  ft  long. 

A channel  will  be  excavated  tlirough  the  saddle  between  the  Yuba  River 
pool  and  the  Dry  Creek  pool  to  facilitate  the  exchange  of  water  between 
the  pools  during  low  lake  levels.  The  lake  will  impound  916,000  acre-ft 
when  the  water  is  at  the  gross  pool,  el  560.  At  this  elevation  the  pool 
will  cover  661tO  acres  and  will  have  67  miles  of  shoreline.  At  minimum 
pool,  el  1(28,  the  lake  will  have  a storage  capacity  of  273,000  acre-ft 
and  a surface  area  of  3l80  acres.  The  afterbay  will  have  a storage 
capacity  of  i*0,i400  acre-ft  and  a surface  area  of  about  9^7  acres  at 
gross  pool,  el  233. 

7.  A 900-megawatt  (MW)  peaking  power  plant  with  four  150-MW  pump 
turbines  and  two  150-MW  conventional  turbines  is  planned  for  the  Yuba 
River  dam.  Water  would  be  released  through  the  main  power  plant  on  an 
intermittent  basis  to  produce  power  during  peak  demand  hours.  When 
power  demand  is  low,  the  pump  turbines  would  pump  water  from  the 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  is  presented  on  page  U. 

**  All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level. 
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Flrure  1.  Location  map 


aftei'bay  to  the  lake  so  the  water  can  be  reused.  Studies  indicate  a 
20-1’W,  base-load  powerplant  downstream  of  the  afterbay  dam  would  be 
feasible. 

8.  Maximum  downstream  temperatures  at  the  mouth  of  the  Yuba  River 
have  been  established  by  the  U.  S.  Fish  and  Wildlife  Service  and  the 
California  Department  of  Fish  and  G.'une  to  maintain  anadromous  fish  runs. 
Minimum  temperatures  have  been  established  for  study  purposes  by  the 

U.  S.  Army  Engineer  District,  Sacramento  (SPK).  Special  intake  struc- 
tures will  be  incorporated  to  enable  withdrawal  of  water  from  different 
levels  in  the  reservoir  to  meet  the  temperature  objectives. 

Approach 

9.  The  thermal  characteristics  of  lakes  are  affected  by  flow 
processes  such  as  inflow  mixing  and  placement,  outflow  withdrawal  dis- 
tribution, diffusion,  and  internal  currents.  For  a lake  having  the 
highly  dynamic  flow  conditions  and  flow  magnitudes  associated  with  the 
pumped-storage  hydropower  activities  of  the  Marysville  Lake  project,  an 
understanding  of  these  flow  processes  is  necessary  if  a reasonable  de- 
scription of  the  thermal  characteristics  is  to  be  made.  A distorted- 
scale  model  (l;l600  horizontal,  l:l60  vertical)  was  used  to  simulate  the 
reservoir  and  afterbay  and  determine  the  response  to  dynamic,  unsteady- 
state,  density-stratified  flow.  Information  from  the  physical  model 
was  used  to  modify  existing  algoritiims  and  develop  new  ones  for  a one- 
dimensional (vertical)  mathematical  model  capable  of  simulating  and 
budgeting  physical  and  conservative  chemical  water  quality 
characteristics . 

10.  The  mathematical  model  allowed  simulation  of  the  hydrodynamic 
and  heat  exchange  characteristics  so  the  thermal  regimes  within  and 
downstream  of  the  project  could  be  evaluated  for  various  hydrologic  and 
meteorologic  conditions. 

11.  Prior  to  the  construction  and  testing  of  the  physical  model, 
preliminary  mathematical  simulations  of  the  Marysville  project  were 
performed  and  a preliminary  report  (App’  ndix  A)  was  furnished  to  GPK. 


The  purpose  of  these  preliminary  simulations  was  to  estimate  the  feasi- 
bility of  the  project  to  meet  downstream  temperature  objectives  and  to 
indicate  port  locations  that  would  enable  these  objectives  to  be  satis- 
fied. The  results  of  the  physical  model  and  refined  mathematical  model 
were  later  used  to  more  accurately  determine  the  ability  of  the  Marys- 
ville project  to  meet  downstream  temperature  objectives. 


PART  II:  PHYSICAL  MODEL 


Purpose 


12.  The  purpose  of  the  distorted-scale  model  was  to  aid  in  defin- 
ing the  hydrodynamics  of  Marysville  Lake  resulting  from  the  planned 
operation  of  the  prototype.  The  model  was  needed  to  determine  the 
effects  of  unsteady  generation  and  pumpback  on  the  density  stratifi- 
cation of  Marysville  Lake  and  afterbay  and  to  help  identify,  validate, 
and  quantify  any  modifications  needed  to  the  mathematical  model  for 
improving  the  reliability  of  the  predictions. 


Scale  Relations 


13-  The  predominant  forces  affecting  density-stratified  flows  in 
lakes  are  inertia  and  gravity  as  modified  by  density  differences.  In 
such  cases,  hydraulic  similarity  between  a model  and  prototype  system 
requires  that  the  ratio  of  inertial  to  gravitational  forces,  defined  as 
the  Froude  number  of  flow,  be  the  same  in  both  the  model  and  the  proto- 
type. With  the  density  differences  in  the  model  set  equal  to  those  in 
the  prototype,  the  accepted  equations  of  hydraulic  similitude,  based  on 
the  Freudian  relations,  were  used  to  express  the  mathematical  relations 
between  the  dimensional  and  hydraulic  quantities  of  the  model  and  the 
prototype.  Allowing  for  vertical  scale  distortion,  the  general  re- 
lations for  transfer  of  model  data  to  prototype  equivalents  are  as  1 


follows : 

Dimension 

Ratio 

Scale  Relation 

Length  in 

vertical  direction 

L = L 

r y 

l:l60 

Length  in 

horizontal  direction 

L = L 
r X 

l:l600 

Area  in  vertical  plane 

A = L L 
r X y 

1:256,000 

Area  in  a 

horizontal  plane 

A = l2 
r X 

1:2,560,000 

(Continued) 
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Dimension 


Ratio 


Velocity 

Time 

Flow 

Density  difference 


V = 

r y 

T = L 
R X y 

Q = L 

r X y 


Ap  = 1 
r 


Scale  Relation 
1:12.65 

1:126.5 

1:3,238,172 

1:1 


Measurements  of  flow,  water-surface  elevations,  and  time  can  be  trans- 
ferred quantitatively  from  the  model  to  the  prototype  by  means  of  the 
scale  relations  above. 


Description 

lU.  The  model  (Figure  2)  was  constructed  to  a distorted  length 
scale  ratio  of  l:l60  vertically  and  l:l600  horizontally.  Through  verti- 
cal scale  distortion,  it  was  possible  to  preserve  turbulent  flow  while 
simulating  the  entire  reservoir  and  afterbay.  It  is  necessary  to  main- 
tain the  same  fundamental  character  of  flow  in  the  model  as  in  the 
prototype.  Use  of  an  undistorted-scale  model  of  the  entire  system  would 
require  a model  of  such  large  dimensions  and  discharge  capacity  to  in- 
sure turbulent  flow  that  the  cost  would  be  Impractical.  By  horizontal 
scale  compression,  the  total  size  of  the  model  was  reduced  and  hydraulic 
similitude  was  preserved. 

15.  The  model  lakes  and  afterbay  were  constructed  of  transparent 
plastic  to  facilitate  photography  and  visual  obse* rations  of  currents 
vertically  and  longitudinally.  The  model  approximated  the  geometry  and 
reproduced  the  scaled  elevation-storage  relationships  of  the  prototype 
lakes.  The  sides  of  the  model  flume  were  stepped  vertically  (Figure  3) 
to  satisfy  the  above-mentioned  requirements  without  visual  distortion 
through  the  side.  Urethane  foam  was  contoured  to  simulate  the  topogra- 
phy upstream  of  the  intake  structures. 

16.  Galine  and  fresh  waters  were  used  to  reproduce  the  density 
variations  that  are  anticipated  in  the  prototype  due  to  temperature 
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Figxire  3.  Example  of  model  design  for  convenient  observation 


I 


differences.  Density  measurements  were  obtained  with  conductivity  and 
temperature  sensors.  The  conductivity  and  density  values  were  converted 
to  density  by  means  outlined  in  Reference  1.  Velocity  distributions 

were  determined  from  video  recordings  of  dye  streak  displacement.  Vari-  ' 

ous  dyes  were  used  to  follow  the  water  movement  of  particular  interest. 

17.  The  preliminary  mathematical  simulations  indicated  that  five  ^ 

vertical  locations  of  selective-withdrawal  ports  el  5^+0,  500,  i4l40, 

375,  and  350)  and  a floodgate  located  near  the  bottom  (invert  el  205) 
could  be  used  to  meet  downstream  temperature  objectives.  The  prototype 
structure  would  consist  of  six  intakes  at  each  vertical  level  to  ac- 
commodate the  six  turbines.  Preliminary  design  established  selective 
withdrawal-port  dimensions  of  50  ft  wide  by  10  ft  high  and  floodgate 
dimensions  of  28  ft  wide  by  28  ft  high.  Because  of  the  difficulty  of 
reproducing  the  small-scaled  horizontal  dimension  of  the  outlet  ports 
(0.031  ft  per  port),  the  horizontal  ports  at  each  level  were  combined 
into  a single  port  in  the  model.  The  width  of  this  port  corresponded 
to  the  width  of  four  ports  which  more  accurately  represented  pumpback 
from  t?ie  four  pump  turbines. 

18.  A width  corresponding  to  six  ports  would  have  better  repre- 
sented flow  during  generation  through  the  six  turbines.  However,  the 
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correct  representation  of  port  width  for  pumpback  was  believed  to  be 
more  critical  than  correct  representation  of  port  width  for  generation. 

A physical  model  test  was  conducted  to  determine  if  the  smaller  width 
would  alter  the  selective  withdrawal  characteristics  during  generation. 
Releases  were  made  through  0.188- ft-  and  0. 125-ft-wide  ports.  There  was 
no  measurable  difference  in  the  outflow  density  of  water  released  from 
the  two  ports  of  different  width.  Therefore,  the  smaller  width  did  not 
alter  the  selective  withdrawal  characteristics.  The  floodgate  was  also 
sized  to  correspond  to  the  width  of  four  ports  and  fixed  in  the  model 
at  invert  el  205  and  centerline  el  219.  The  five  selective  withdrawal 
port  levels  were  simulated  with  a single  movable  port,  which  represented 
four  intakes  at  a given  level,  that  could  be  positioned  at  any  one  of 
the  five  port  elevations.  This  upper  movable  port  emptied  into  a wet- 
well  and  was  positioned  by  an  electro-mechanical  actuator. 

19*  Generation  and  pumpback  hydrographs  were  simulated  using  a 
variable  speed,  reversible,  programmable,  positive  displacement  pump. 
Actuated  butterfly  valves  were  used  to  direct  the  flow  through  either 
the  upper  (selective  withdrawal)  port  or  the  lower  (floodgate)  port.  A 
constant  head  tank,  rotameter,  and  gate  valve  were  used  to  provide  and 
regulate  inflow  from  the  Yuba  River.  Downstream  releases  from  the  after- 
bay were  controlled  by  a rotameter  and  gate  valve. 

20.  Model  control  devices  could  be  operated  manually  or  automati- 
cally from  the  control  console  (Figure  2b).  Automatic  control  was 
achieved  through  analog  signals  generated  by  magnetic  card  readers  that 
tracked  input  data.  Automatic  control  allowed  unsteady  operation  to 
proceed  for  extended  simulation  periods  (It  weeks  prototype  time).  Be- 
cause the  model  was  not  able  to  reproduce  surface  heat  exchange,  real- 
istic prototype  simulations  exceeding  ^t  weeks  were  not  attempted  because 
the  meteorological  effects  could  exceed  the  hydrodynamic  effects. 

Model  Tests  and  Results 


Prototype  simulation 

21.  The  anticipated  thermal  stratification  conditions 
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(see  Appendix  A)  for  July  of  an  average  hydrologic  and  meteorologic  year 
(1962)  were  simulated  in  the  model.  The  model  simulated  a constant  in- 
flow from  the  Yuba  River  of  2288  cfs  at  a temperature  of  lU.7°C.  The 
planned  project  operation  scheme  for  the  above-mentioned  average  year, 
which  was  furnished  by  the  Sacramento  District  and  was  applied  in  the 
model,  is  presented  below: 


Turbine  Configuration 

2 conventional  and 
reversible  turbines 


Generation  Flow 

37,700  cfs  for 
3.71  hr  per  day, 
6 days  per  week; 
start  at  5 pm 


Afterbay 

Pumpback  Flow  Release 

8200  cfs  for  2360  cfs 

9 hr  per  night,  constant 

6 days  per  week;  release 

start  at  10  pm 


Generation  and  pumpback  occur  Monday  through  Saturday.  With  the  pool  at 
el  53^,  the  upper  level  port  center  line  was  set  at  el  500  for  with- 
drawal during  model  simulation.  Withdrawal  through  this  port  would  be 
required  to  meet  downstream  temperature  objectives  during  July  of  an 
average  year.  Two  pumpback  configurations  were  tested,  pumpback  through 
the  withdrawal  port  and  pumpback  through  the  low  level  flood  gate. 
Marysville  Lake  stratification 

22.  The  effects  of  unsteady  operation  on  Marysville  Lake  stratifi- 
cation were  investigated  by  operating  the  hydrodynamic  model  with  the 
anticipated  prototype  operation  scheme.  Changes  in  the  stratification 
were  evaluated  by  comparing  vertical  density  profiles  taken  at  the  end 
of  each  prototype  week  of  simulation.  As  anticipated,  these  density 
profiles  revealed  that  the  pumpback  flow  process  had  a much  stronger 
effect  on  density  stratification  than  did  inflow  or  withdrawal.  Pumped 
flow  establishes  a density  current  that  causes  mixing  with  respect  to 
time  within  a vertical  zone  of  the  pool  as  shown  in  Figure  U.  This  same 
phenomenon  was  observed  in  a previous  pumped-storage  model  study  and  is 
discussed  in  detail  in  Reference  2.  This  zone  of  mixing  corresponded 

to  the  zone  of  interflow  of  the  pumpback  current.  The  observed  puinpback 
flow  process  is  schematically  described  in  F’igure  5. 

23.  To  account  for  the  mixing  within  the  zone  of  interflow  of  the 
pumpback  current,  a procedure  that  was  developed  during  the  previous 
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Figure  5.  Puinpback  characteristics 
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study  was  found  to  be  appropriate  for  describing  the  effect  in  Marys- 
ville Lake.  The  mixing  procedure  consists  of  three  parts:  prediction 

of  the  pumpback-current  thickness,  D,  in  the  vicinity  of  the  dam;  de- 
termination of  mixing  coefficients;  and  partial  mixing  of  layers  within 

the  zone  of  the  current.  Mixing  coefficients  characteristic  of  the  i 

Marysville  project  were  determined  from  density  profiles  taken  in  the  j 

model  with  respect  to  time.  The  mixing  technique  is  discussed  in  i 

greater  detail  in  Part  III. 

Marysville  generation  releases 

2h.  The  densities  of  flows  released  into  the  afterbay  were  moni- 
tored during  simulated  power-generation  operations  to  evaluate  whether 
variations  of  release  density  resulted  from  the  unsteady  operation. 

There  were  no  transient  effects  on  release  densities  observed  during 
any  one  generation  period.  Release  densities  varied  some  from  one 
period  to  the  next  but  not  during  a period.  The  observed  release  densi- 
ties compared  closely  with  those  computed  by  the  U.  S.  Army  Engineer 

Waterways  Experiment  Station  (WES)  generalized  selective  withdrawal 
3 

technique.  It  was  necessary  to  use  the  actual  discharge  rate  that  i 

occurred  within  a generation  period  rather  than  a daily  average  gener-  i 

ation  rate  to  make  reasonable  predictions.  Actual  generation  rates  are 
used  in  the  mathematical  model  (paragraph  U6) . 

Afterbay  mixing 

25.  Tests  were  conducted  to  determine  the  mixing  characteristics 
of  the  afterbay.  The  afterbay  was  set  up  with  a density  stratification 
corresponding  to  23°C  and  l8°C  waters  for  the  epilimnion  and  hypolimnion, 
respectively,  and  the  operation  schedule  was  imposed.  The  vertical 
density  stratification  of  the  afterbay  did  not  exist  after  one  week  of 

prototype  simulation.  Therefore,  stratification  of  the  afterbay  is  not  ■ 

expected  to  exist.  The  afterbay  was  found  to  be  almost  totally  mixed  at 
the  end  of  each  prototype  day.  However,  it  was  observed  that  during,  the 
period  of  a day  various  stages  of  mixing  existed.  During  generation,  the 
release  waters  from  Marysville  Lake  partially  mixed  with,  but  primarily 
displaced,  the  pregeneration  afterbay  water,  thereby,  occupying  tfie  up- 
stream reach  of  the  afterbay.  Because  there  was  some  mixing,  the 
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pumpback  water  consisted  of  pregeneration  afterbay  water  and  Marysville 
Lake  release  water.  After  the  completion  of  pumpback,  the  remaining 
generation  water  was  observed  to  disperse  throughout  the  afterbay.  For 
numerically  modeling  the  afterbay,  this  allowed  the  simplifying  assump- 
tion that  the  afterbay  is  fully  mixed  at  the  end  of  each  day. 

Pumpback  temperature 

26.  Tests  were  conducted  to  determine  what  volume  of  pregeneration 
afterbay  water  was  pumped  back  with  the  Marysville  release  water  follow- 
ing a generation  period.  This  information  was  needed  to  provide  proper 
simulation  of  the  pumpback  temperature  and  mixing  of  the  residual  waters 
in  the  afterbay  after  pumpback.  Fluorescent  dye,  which  was  mixed  with 
the  pregeneration  afterbay  water,  was  used  to  trace  the  contribution  of 
pregeneration  afterbay  water  in  the  pumpback  flow. 

27.  The  predominant  factor  influencing  the  pumpback  constituents 
was  considered  to  be  the  pumped  volume  ratio,  which  is  defined  as: 


FTunped  volume  ratio  = 


volume  of  water  pumped  back  in  one  day 
volume  of  water  released  during 
generation  in  one  day 


(1) 


To  examine  a full  range  of  flow  conditions,  three  operation  schedules 
were  selected  for  model  testing  from  the  various  operations  anticipated 
by  SPK.  The  flows  and  pumped  volume  ratios  for  the  conditions  tested 


are  presented  below. 

Pumped 

Volume 

Generation  Flow Pumpback  Flow Ratio 

36,  Loo  cfs  for  1.80  hr  per  day  1000  cfs  for  2.h2  hr  per  day  O.OL 

37,700  cfs  for  3-71  hr  per  day  8200  cfs  for  9 hr  per  day  0.53 

37,  Loo  cfs  foi  3.28  hr  per  day  9280  cfs  for  9 hr  per  day  0.68 


A constant  afterbay  release  of  2360  cfs  was  used  for  each  condition.  A 
measurement  of  initial  fluorescent  dye  concentration  in  the  afterbay  was 
made  prior  to  generation.  The  generation  period  was  then  simulated 
causing  mixing  in  the  afterbay.  After  generation,  pumpback  was 


•i 

( 
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initiated.  During  purapback,  samples  of  water  drawn  into  the  draft  tubes 
were  collected  and  analyzed  for  fluorescent  dye  concentration.  The 
average  dye  concentration  of  the  pumped  flow  was  related  to  the  initial 
concentration  of  dye  in  the  afterbay  to  determine  the  volume  of  pregener- 
ation afterbay  water  pumped  back.  The  ratio  of  the  volume  of  pregener- 
ation afterbay  water  pumped  back  with  respect  to  the  total  volume  of 
water  pumped  back  in  one  day  was  correlated  with  the  pumped  volume  ratio. 
A best  fit  of  the  data  gave  the  following  equation: 

V = V (0.00i»e^'®°^)  (2) 

A p 


where 


V 

P 

e 


r 


volume  of  pregeneration  afterbay  water  that  is  piimped  back  in 
one  day,  acre- ft 

total  volume  of  water  that  is  pumped  back  in  one  day,  acre-ft 
natural  logaritiimic  base,  2.7l8j 
pumped  volume  ratio 


With  a knowledge  of  V^,  it  is  possible  in  the  numerical  model  to  detei-- 
mine  the  pumpback  temperature  and  mix  the  water  remaining  in  tiie  after- 
bay after  pumpback. 

Pumpback  entrainment 

28.  The  amount  of  flow  entrained  by  the  pumpback  jet  was  deter- 
mined from  the  piiysical  model  for  both  the  high-  and  low-level  pumping 
configurations.  Entrained  flow  is  water  that  is  j .lied  into  and  mixed 
with  the  pumpback  Jet.  The  mixture  goes  into  storage  in  the  pool  ;it  a 
density  equal  to  the  density  of  the  mixture,  thus,  creating  a pumpback 
density  current.  Coefficients  of  entrainment  are  necessary  for  numericai 
model  simulations. 

29.  Conservation  of  mass  and  volume  were  used  to  deterriine  t(ie 
amount  of  bulk  entrainment  as  follows.  From  the  conservation  of  ma.m 


where 


p = averat'e  density  of  pum])back  current,  g/cc 
= volume  flow  rate  of  pumpback  current,  cfs 
= density  of  pumpback  jet,  g/cc 
= volume  flow  rate  of  pumpback  Jet,  cfs 

p = average  density  of  entrained  current,  g/cc 
e 

Q = volume  flow  rate  of  entrained  current,  cfs 
e 


Also,  from  the  continuity  of  volume  for  an  incompressible  fluid 


Q = Q + Q 
c o e 


Substitution  of  equation  1 into  equation  3 for  yields 


Q = EQ 
e o 


where  E is  the  entrainment  coefficient  and 


P - 
o c 


c e 


The  value  of  p^  was  taken  from  the  initial  density  profile  at  the  ele- 
vation where  the  maximum  contribution  of  entrainment  was  observed  in  the 
model.  Density  measurements  were  taken  at  the  Marysville  tailrace  and 

within  tlie  pumpback  current  to  olitain  p and  p , respectively. 

o c 

30.  With  pumpback  through  an  upper  port,  E was  found  to  be  approxi- 
mately 0.7  for  the  conditions  tested.  The  entrainment  was  observed  to 
occur  primarily  from  ttie  layers  of  the  pool  at  and  immediately  above  tiie 
location  of  the  port.  Considerably  more  entrainment  was  detected  witti 
pumpback  through  the  lower  level  floodgates.  Calculations  from  density 
measurements  indicated  that  E was  as  large  as  2.5  to  3.0.  To  substanti- 
ate tiiese  values,  velocity  profiles  of  tlie  entrained  flow  were  measured. 
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Assuming  uniform  flow  laterally,  the  velocity  profiles  were  converted  to 
flow  rate,  or  Q . By  substitution  of  Q and  Q into  equation  5,  E was 

G GO 

found  to  be  about  2.5.  For  pumpback  through  the  floodgate,  entrainment 
occurred  from  the  bottom  up  to  an  elevation  about  170  ft  above  the 
bottom. 

Inflow 

31.  With  a constant  inflow  from  the  Yuba  River  of  2288  cfs  simu- 
lated in  the  model,  the  inflow  current  traveled  to  the  Marysville  Dam 

in  a week  (prototype  time).  This  travel  time  indicates  that  the  thermal 
structure  of  the  lake  at  the  dam  should  respond  in  a short  time  to  the 
effects  of  inflow.  Practically  no  entrainment  due  to  inflow  was  ob- 
served in  the  model  for  the  2288-cfs  flow  rate.  Therefore,  the  assump- 
tion of  zero  inflow  entrainment  was  used  for  the  numerical  simulations. 
Topographical  effects 

32.  The  ridge  upstream  of  and  generally  parallel  to  the  dam  was 
observed  to  influence  circulation  and  mixing  in  the  vicinity  of  the  dam 
during  pumpback  through  the  floodgate.  For  both  pumpback  configurations 
(jet  entering  the  pool  through  an  upper  level  port  or  through  the  low 
level  floodgate),  the  pumpback  current  travelled  upstream  and  spread 
throughout  Marysville  Lake.  For  pumpback  through  an  upper  level  port 
with  the  conditions  simulated  in  the  physical  model  (paragraph  21 ) , the 
current  was  able  to  spread  over  the  top  of  the  ridge  as  it  travelled 
upstream.  For  pumpback  through  the  floodgate,  the  plume  became  neutrally 
buoyant  at  an  elevation  within  the  water  column  (el  1+90-500)  that  was 
below  many  of  the  ridge  peaks.  This  forced  much  of  the  current  to 
travel  around  Timbuctoo  Bend  rather  than  over  the  ridge.  Additionally, 
the  entrained  flow,  which  consisted  of  water  between  the  bottom  and 
about  el  370,  was  forced  to  travel  around  Timbuctoo  Bend  because  of  tlie 
presence  of  the  ridge.  For  the  low  level  pumpback,  the  ridge  certainly 
influenced  circulation  in  the  vicinity  of  the  dam,  and  the  flow  restric- 
tion created  by  the  ridge  tended  to  increase  the  turbulence  and  thick- 
ness of  the  pumpback  current  in  the  region  between  the  d;un  and  the 
ridge.  The  above-mentioned  phenomena  in  addition  to  the  large  entrained 
flow  that  occurs  with  the  low  level  pumpiiack  of  warm  water  would  cause 


22 


some  mixing  and  warming  of  water  in  the  hypolimnlon.  If  an  objective 
of  the  project  is  to  conserve  cold  water  reserves  for  fall  releases, 
then  pumpback  through  the  floodgate  (accompanying  releases  through  an 
upper  port)  could  he  undesirable  and  pumpback  at  the  level  of  withdrawal 
is  recommended. 

33.  The  invert  of  the  channel  (el  ^400)  connecting  the  Dry  Creek 
embayment  with  the  main  lake  prevents  water  in  the  Dry  Creek  embaynient 
below  el  400  from  entering  the  main  lake.  This  reduces  the  total  cold 
water  storage  capacity  of  the  lake.  To  account  for  this  reduction  in 
the  numerical  model,  the  voluines  of  layers  below  el  hOO  that  are  con- 
tained in  the  Dry  Creek  embayment  were  subtracted  from  the  total  volume 
of  each  layer  below  el  UOO.  Although  water  below  el  UOO  in  the  Dry 
Creek  embayment  was  stagnant,  water  above  el  1+00  circulated  freely. 
Inflow,  pumpback,  and  witndrawal  flows  passed  to  and  from  the  Dry  Creek 
embayment  above  el  1+00. 
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PART  III:  MATHEMATICAL  MODEL  DESCRIPTION 


3^ . The  internal  heat  budget  in  Marysville  Lake  and  the  release 
temperatures  downstream  of  the  project  and  at  the  confluence  of  the  Yuba 
and  Feather  Rivers  were  predicted  using  a numerical  simulation  model. 

The  model  used  in  conjunction  with  this  investigation  was  previously 
developed  at  WES  and  modified  to  uniquely  describe  the  thermal  character- 
istics of  the  Marysville  Lake  project.  During  the  development  of  the 
model,  considerable  insight  has  been  gained  from  the  research  efforts  of 
Bohan  and  Crace,  Clay  and  Fruh,  Edinger  and  Geyer,  Dake  and  Harleman, 
and  others. 


Fundamental  Assumptions 


35-  The  Marysville  Lake  hydrothermal  mathematical  model  provides 
a procedure  for  examining  the  balance  of  thermal  energy  imposed  on  an 
impoundment.  This  energy  balance  and  lake  hydrodynamic  phenomena  are 
used  to  map  vertical  profiles  of  temperature  and  release  temperatures  in 
the  time  domain.  The  model  includes  computational  methods  for  simulat- 
ing heat  transfer  at  the  air-water  interface,  advective  heat  due  to  in- 
flow, outflow,  and  pumpback  processes,  and  the  internal  dispersion  of 
thermal  energy.  The  model  is  conceptually  one-dimensional  based  on  the 
division  of  the  impoundment  into  discrete  horizontal  layers  of  uniform 
thickness.  Assumptions  include  the  following: 

a.  Isotherms  are  parallel  to  the  water  surface  both  laterally 
and  longitudinally. 

b.  The  water  in  each  discrete  layer  is  isotropic  and  physi- 
cally homogeneous . 

c_.  Internal  advection  and  heat  transfer  occur  only  in  tiie 
vertical  direction. 

External  advection  (inflow,  outflow,  and  pumpback)  occurs 
as  a uniform  distribution  within  each  layer. 

£.  Internal  dispersion  (between  layers)  of  thermal  ener/iy  is 
accomjilished  by  a diffusion  mechanism  that  combines  the 
effects  of  molecular  diffusion,  turbulent  diffusion,  and 
thermal  convection. 


36.  The  surface  heat  exchange,  internal  mixin^^,  and  advection 
processes  are  simulated  separately  and  their  effects  are  introduced 
sequentially  at  specified  time  intervals.  A simplified  flow  chart  of 
the  mathematical  simulation  procedure  is  presented  in  Fip-ure  6. 

Surface  Heat  Exchanpie 

37.  The  Maj-ysville  L<ike  .mathematical  model  employs  an  approach  to 
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the  evaluation  of  surface  heat  transfer  developed  by  Edinger  and  Geyer. 
This  method  formulates  equilibrium  temperatures  and  coefficients  of 
surface  heat  exchanre.  Equilibrium  temperature  is  defined  as  that 
temperature  at  which  the  net  rate  of  heat  exchange  between  the  water 
surface  and  the  atmosphere  is  zero.  The  coefficient  of  surface  heat 
exchange  is  the  rate  at  which  the  heat  transfer  process  will  occur.  The 
equation  describing  this  relationship  is: 

H = K(E  - 8)  (T) 


where 

2 

H = net  rate  of  surface  heat  transfer,  Btu/ft  /day 

? o 

K = coefficient  of  surface  heat  exchange,  Btu/ft  /day/  F 
E = equilibrium  temperature,  °F 
8 = surface  temperature,  °F 

The  computation  of  equilibrium  temperature  and  heat  exchange  coefficient 

is  ijased  solely  on  meteorological  data  as  outlined  by  Edinger, 
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Duttweiler,  and  Geyer. 

38.  Tiie  net  iieat  exchange  at  the  surface  is  composed  of  seven 
heat  exchange  processes: 

a_.  Shortwave  solar  radiation. 

Reflected  shortwave  radiation. 
e_.  Long-wave  atmospiieric  radiation, 
d.  Reflected  long-wave  radiation. 

£.  Heat  transfer  ilue  to  conduction, 
f.  Back  radiation  from  the  water  surface. 


Heat  loss  due  to  evaporation. 


For  every  day  of  meteorological  data,  the  seven  heat  exchange  terms  can  < 

be  evaluated  and  the  net  heat  exchange  expressed  in  terms  of  an  equilib- 
rium temperature  and  an  exchange  coefficient. 

39.  All  of  the  surface  heat  exchange  processes,  with  the  exception 
of  shortwave  radiation,  affect  only  approximately  the  top  few  feet  of 
the  lake.  Shortwave  radiation  penetrates  the  water  surface  and  increases 
the  temperature  at  greater  depths.  Based  on  laboratory  investigations, 

Dake  and  Harleman^  have  suggested  an  exponential  decay  with  depth  for 
describing  the  heat  flux  due  to  shortwave  penetration. 

ho.  The  surface  heat  exchange  concepts  are  implemented  in  the 
model  by  the  exponential  penetration  of  a percentage  of  the  incoming 
shortwave  radiation  and  the  placement  of  the  effect  of  all  other  sources 
of  surface  heat  exchange  into  the  surface  layer.  This  procedure  can  be 
expressed  mathematically  by  the  following  two  equations: 

H = K(E  - Q)  - (1  - 3)S  (8) 

s 

H.  = (1  - 6)Se"^*i  (9) 

where 

2 

= heat  transfer  rate  into  or  out  of  surface  layer,  Btu/ft  /day 

6 = shortwave  radiation  absorbed  in  the  surface  layer,  percent 

2 

S = total  incoming  shortwave  radiation  rate,  Btu/ft  /day 

2 

= rate  of  heat  absorbed  in  layer  i , Btu/ft  /day 
e = natural  logarithmic  base  (2.7183) 

A = heat  absorption  coefficient,  ft  ^ 

= depth  below  surface,  ft 

Ul.  Equations  8 and  9 are  applied  to  the  Marysville  Lake  temper- 
ature profile  once  during  each  one-day  time  step.  The  net  heat  exchange 
rate  into  each  layer  is  computed  and  converted  to  a temperature  change. 

The  temperature  changes  are  used  to  determine  an  uf'dated  temperature 
profile  for  Marysville  Lake.  Numerically,  the  Marysville  afterbay  is 
mixed  at  the  end  of  each  daily  time  step.  Thu.s , surface  heat  exchange 
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in  the  afterbay  is  effected  by  converting  the  total  heat  exchange  rate 
as  indicated  in  equation  T to  a temperature  change.  The  temperature 
change  is  then  added  to  the  existing  afterbay  temperature  once  during 
each  one-day  time  step. 


Inflow 


i+2.  The  process  of  inflow  into  a lake  is  simulated  in  the  mathe- 
matical model  by  the  placement  of  inflow  quantity  and  quality  at  tiiat 
layer  in  which  the  density  of  the  lake  corresponds  most  nearly  to  the 
density  of  the  inflow.  Although  not  indicated  in  this  study,  research 
efforts  and  physical  model  studies  at  WES  have  shown  that  entrainment- 
induced  density  currents  can  exist  which  flow  upstrear;  along  the  surface 
into  the  turbulent  mixing  zone  caused  by  the  inflow. 

^+3.  When  it  exists,  entrainment  is  implemented  in  t)ie  model  by 
augmenting  the  inflow  quantity  with  a volume  from  the  surface  layer. 

The  volume  from  the  surface  layer  is  generally  expressed  as  a percentage 
of  the  inflow  quantity.  Characteristics  of  inflow  and  the  entrained 
flow  are  averaged,  and  mixed  values  of  density,  temperatuia  , an<l  other 
water-quality  parameters  are  determined.  The  mixed  density  is  useii  to 
determine  placement  of  the  total  quantity  and  mixed  quality.  Gimulation 
of  the  inflow  process  displaces  upward  a volume  equal  to  the  total  in- 
flow quantity.  This  upward  displacement  is  reflected  in  the  model  ly 
an  increase  in  the  water  surface.  A corresponding  decrease  in  water 
surface  occurs  as  a result  of  the  outflow  simulation  process. 

Internal  Dispersion 

hh.  The  internal  dispersion  process  is  represented  by  an  internaJ 
mixing  scheme  based  on  a simple  diffusion  analogy.  Internal  mixing 
transfers  heat  and  other  water-quality  constituents  between  ad.iaeent 
layers.  The  magnitude  of  the  transfer  between  two  layers  is.  a perceii' - 
age  of  the  total  transfer  required  to  completely  mix  tiie  two  layers. 

This  percentage  is  a mixing  coefficient  wiilcli  is  d('fined  for  every 
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layer.  Data  input  includes  values  of  the  mixing  coefficient  at  the  top 
and  at  the  bottom  of  the  lake.  An  exponential  fit  between  the  two  ex- 
treme values  is  used  to  determine  the  appropriate  coefficient  for  each 
layer . 

Outflow 

I45.  Tlie  outflow  component  of  the  model  incorporates  the  selective- 

3 

Withdrawal  techniques  for  orifice  flow  developed  at  WES.  Additionally, 
the  selective-withdrawal  description  was  analyzed  to  include  the  appro- 
priate application  to  simultaneous  releases  from  multiple  horizontal 
ports  (Appendix  B).  Transcendental  equations  defining  the  location  of 
zero- velocity  limits  are  solved  with  a half-interval  search  method.  The 
location  of  the  zero-velocity  limits  is  functionally  dependent  on  the 
coiifiguration  of  the  withdrawal  device,  release  flow  rate,  and  density 
structure  within  the  lake.  With  knowledge  of  the  limits  of  the  with- 
drawal zone,  the  velocity  profile  due  to  outflow  can  be  determined.  The 
flow  from  each  layer  is  then  computed  as  the  product  of  the  velocity  in 
the  layer  and  the  width  and  thickness  of  the  layer.  A flow-weighted 
average  is  applied  to  the  temperature  profile  to  determine  the  value  of 
tiie  release  temperature  for  eacii  flow  condition.  When  the  port  to  be 
used  for  outflow  is  not  specified,  such  as  in  a predictive  mode,  the 
model  compares  target  terap'eratures  and  the  simulated  tiiermai  profiles 
to  select  the  ports  to  be  used  for  withdrawal. 


Operation  Scheduler. 


1*6.  liecause  tiie  generation  and  pumpback  flows  cliange  markedly  in 
a day  and  the  rate  of  these  flows  affects  witlidrawal  and  pumpback 
ciiaracteri sties , it  was  not  adequate  to  use  daily  average  flows  to  accu- 
rately simulate  generation  and  pumpback  op'crations.  Daily  average  flow 
routings  usually  are  quite  ade(iuate  for  less  dynamic  reservoirs,  such 
as  water  supjily  and  flood  control  pro.lect:;.  For  this  model,  tiie  I'lanneii 
generation  and  pumpback  flow  rates  and  dui’at,  ions  (in  liours)  are  us('d 
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rather  than  daily  averages.  The  model  can  handle  up  to  twelve  difl'erent 
flow  conditions  within  a day.  The  .Tulian  day  number,  t)ie  modes  of  oper- 
ation (generation  or  pumpback),  flow  rates,  durations  of  flow,  and  ele- 
vation of  tlie  pumpback  device  are  input  for  ■ ach  day  that  a change  in 
the  operation  scliedule  occurs. 

Marysville  Afterbay 

^7.  The  afterbay  is  modeled  numerically  by  maintaining  heat  and 
water  budgets.  The  afterbay  is  fully  mixed  at  the  end  of  each  daily 
time  step.  Tliis  assumption  is  supported  by  physical  model  observations 

g 

and  densimetric  Froude  number  calculations,  and  it  allows  further 
simplifying  assumptions  to  be  made  with  respect  to  advection  and  lieat 
exchange  in  the  afterbay.  The  daily  surface  heat  exchange  is  accom- 
plished by  applying  equation  T to  the  afterbay  surface  area  that  exists 
at  the  beginning  of  each  daily  time  step.  Once  during  each  day  of  simu- 
lation, the  volume  of  water  in  the  lake  is  adjusted  to  sequentially 
account  for  (a)  Marysville  Lake  generation  volume  that  is  not  pumped 
back,  (b)  pregeneration  afterbay  volume  that  is  pumped  back,  and  (c) 
afterbay  volume  tliat  is  released  downstream.  The  net  contributions  of 
volume  from  Marysville  Lake  at  the  res]'cctivo  temperatures  are  added  to 
the  afterbay  volume  producing  a mixed  value  of  temperature.  This  mixed 
temperature  is  the  afterbay  release  temperature  for  that  day  of  simu- 
lation. Figure  7 helps  describe  the  approacii  followed  in  modeling  tiie 
afterbay . 

Pumpback  Temperature  and  Kntrainnient 


U8.  For  the  numerical  simulations,  it  is  necessary  to  define  tiie 
temperature  of  tlie  flow  punqied  from  the  afterbay  to  tlie  lake.  Tf  gener- 
ation does  not  occur  on  the  same  day  that  tiiere  is  pumpback,  then  ttic 
temperature  of  tlie  afterbay  water  constitutes  the  pimiped  flow  ti'mper- 
ature.  However,  if  generation  occurs  prior  to  pumpback  on  the  same  day, 
then  the  pumped  temperature  can  result  from  a mixture  of  pregeneration 
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afterbay  water  and  water  released  during  generation.  Equation  2 is  used 
to  determine  what  amounts  of  the  pumped  water  are  pregeneration  afterbay 
and  previously  generated  water.  The  temperature  of  the  pumped  flow  is 
determined  from  a volume-weighted  average  of  the  constituents.  The 
residual  water  of  the  afterbay  after  pumpback  is  mixed  to  obtain  a 
volume- weighted  average  temperature. 

U9.  Entrainment  is  expressed  as  a percentage  of  pumpback  water 
volume.  The  entrained  water  of  the  Marysville  pool  at  its  temperature 
is  withdrawn  and  numerically  mixed  with  the  pumpback  volume  at  its 
temperature.  This  total  volume  of  water  at  the  mixed  temperature  is 
placed  into  the  layer  of  the  lake  that  most  closely  corresponds  to  the 
mixed  temperature  of  pumpback.  For  pumpback  through  the  upper  ports,  an 
entrainment  coefficient,  E,  of  O.Y  was  determined  from  the  physical 
model  tests  and  was  used  in  the  numerical  simulation  model.  The  layer 
of  the  pool  from  which  entrainment  occurs  is  defined  as  the  layer  immedi- 
ately above  the  layer  of  the  pumpback  port.  For  pumpback  through  the 
floodgate,  a value  of  2.5  was  used  for  E (see  paragraph  30)  and  water 
was  entrained  from  the  elevation  of  the  center  line  of  the  floodgate. 

Pumpback  Mixing  Technique 

50.  From  the  piiysical  model  observations  it  was  noted  that  differ- 
ent density  profiles  existed  in  the  vicinity  of  tlie  d.im  relative  to 
those  observed  furtiier  upstream  in  the  lake.  Because  the  mathematical 
model  is  one-dimensional,  the  decision  was  made  to  simulate  the  tliermal 
profile  immediately  upstream  of  the  dam,  since  this  profile  controls  tfie 
withdrawal  and  pumpback  characteristics  within  the  lake. 

51.  In  the  matiiematical  model,  the  pumpback  current  is  placeii  in 
the  layer  whose  density  most  nearly  corresponds  to  tiie  density  of  the 
pumpback  current.  In  the  vicinity  of  the  dam,  the  pujnpluick  current 
actually  causes  mixing  in  a r.one  about  tliis  layer  ar,  inciicated  by  the 
physical  model.  To  rej'rttsent  this  mixing  effect,  and  still  maintain  Uu' 
heat  budget  of  Uie  lake,  Uie  layi'rs  wiUiin  Uie  r.one  are  numt'rieally 
mixed.  The  mixing  technique  employed  was  develepi>d  for  a p^evi^■u;■. 


pumped-storage  model  study.  The  teclinique  is  based  on  the  concept  of  a 
portion  of  each  of  the  layers  in  the  pumpback  zone  being  removed,  mixed 
together,  returning  tfie  mixed  water  to  the  layers  in  the  zone,  and  then 
mixing  within  each  layer.  The  portion  of  water  removed  and  mixed  from 
each  layer  is  computed  by  multiplying  a mixing  coefficient  (a  number 
between  0 and  l)  times  the  volume  of  the  layer.  The  mixing  coefficients 
were  determined  using  data  from  the  physical  model.  Density  profiles 
were  measured  in  the  physical  model  before  and  after  pumpback  operations. 
By  analyzing  the  difference  in  these  profiles,  it  was  possible  to  deter- 
mine the  mixing  coefficients  necessary  for  application  to  the  initial 
profile  so  that  the  final  profile  could  be  obtained. 

52.  The  thickness,  D (ft),  of  the  pumpback  cm-rent,  which  is  the 
same  as  the  above-mentioned  mixing  zone,  is  calculated  from  the  equation 


where 

Q = volumetric  flow  rate  of  the  pumpback  current  (includes  en- 
trainment),  cfs 

W = average  reservoir  width  at  the  elevation  of  the  current  in 
the  vicinity  of  the  dam,  ft 

Ap  = density  difference  of  the  epilimnion  and  hypolimnion,  g/cc 

p = average  density  of  the  pumpback  current,  g/cc 
c 

2 

This  equation  was  developed  from  the  previous  model  study. 

53.  Mixing  coefficients  were  computed  for  the  pumpback  operation 

tests  conducted  in  the  physical  model  (paragraph  2l).  Analysis  of  the 

mixing  coefficients  revealed  that  they  had  characteristics  similar  to 

> 

the  mixing  coefficients  developed  in  the  previous  model  study. based 
on  those  results,  plus  the  results  from  the  Marysville  Lake  pumpback 
operation  tests,  an  equation  was  developed  to  predict  the  mixing,  coel’- 
ficient  in  a layer  as  a function  of  the  distance  of  that  layer  from 
tfie  pumpback  current  inflow  layer 
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(11) 


"i 


A -BX- 

= A e 1 


where 

= mixing  coefficient 
A = 0.16T5 
B = 0.118 

= distance  from  the  pumpback  current  inflow  layer  to  layer  i 
i = layer 


5I.  To  represent  the  pumpback  mixing  technique  within  the  mathe- 
matical model,  the  following  sequence  of  computational  steps  is  employed. 
The  thickness  of  the  zone  affected  by  the  pumpback  current  is  computed 
and  this  zone  is  distributed  about  the  layer  of  pimipback  placement. 

The  mixing  coefficient  for  each  layer  in  the  zone  is  computed  by  equa- 
tion 11,  the  amount  of  each  layer  to  be  mixed  is  computed  as  the  product 
of  the  mixing  coefficient  and  the  layer  volume,  and  the  portions  from 
the  layers  are  mixed  together.  The  mixed  water  is  returned  to  the 
layers,  and  each  layer  is  individually  mixed. 

Downstream  Temperature  Routing 

55.  Afterbay  temperatures  were  routed  downstream  using  Edinger's^ 
temperature  model  for  a well  mixed  stream  with  uniform  depth  and  ve- 
locity : 

(-KX/yc  vd) 

T = (T  - E)e  ^ + E (12) 

X o 

where 

T^  = temperature  at  downstream  end  of  reacli, 

T^  = temperature  at  upstream  end  of  reach, 

E = equilibrium  temperature, 

2 o 

K = surface  heat  excfiange  coefficient,  Btu/ft  /day/  F 


’’M 
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X = reach  length,  ft 

Y = specific  weight  of  water,  62. h Ib/ft^ 

c = specific  heat  of  water,  1 Btu/lb/°F 

P 

d = uniform  stream  depth,  ft 

V = uniform  stream  velocity,  ft/day 

56.  Routings  were  accomplished  in  two  steps  due  to  projected  irri- 
gation withdrawals  at  Daguerre  Point.  Afterbay  temperatures  were  routed 
i-3  miles  to  Daguerre  Point  with  a velocity  and  depth  corresponding  to 
the  afterbay  release  flow.  The  temperature  computed  for  Daguerre  Point 
was  then  routed  9-2  miles  to  the  mouth  of  the  Yuba  River  with  a velocity 
and  depth  corresponding  to  the  flow  remaining  after  the  diversion. 
Equilibrium  temperatures  and  surface  heat  exchange  coefficients  were 
the  same  as  those  used  for  Marysville  Lake  and  afterbay.  Velocity-depth 
products  were  computed  from  the  relationship 

vd  = (13) 


where 

2 

vd  = velocity  times  depth,  ft  /day 
Q = flow  rate,  cfs 

This  relationship  was  developed  by  SPK  from  fitting  observed  flow  data 
for  water  year  1975  with  values  of  vd  computed  from  equation  12  using 
observed  water  temperature  and  meteorological  data  for  1975- 


PART  TV:  NTATHErTATICAL  SIMULATIONS 


Development  of  Data 

57-  At  the  request  of  SPK,  thj’ee  hydrolo/^ic  years  were  used  for 
the  simulations;  193^<,  the  driest  year  in  the  period  of  record;  19^2, 
the  wettest  year  in  the  period  of  record;  and  1962,  an  average  year  for 
the  period  of  record.  Meteorological  data  from  Beale  Air  Force  Base, 
located  about  7 miles  south  of  the  project,  exist  from  I960  to  the 
present.  Equilibrium  temperatures  were  computed  from  these  meteorologi- 
cal data  for  years  i960  through  1970.  These  equilibrium  temperatures 
were  analyzed  to  determine  hot,  1967;  cold,  1963;  and  average,  1962, 
meteorological  years.  Tliese  hydrologic  and  meteorologic  conditions  were 
combined  to  form  three  study  years  composed  as  follows : 


Study 

Year 

1 

1962 

1962 

Hydrology  (average) 
Keteorolo(-’y  (average 

Study 

Year 

2 

193^* 

1967 

Hydrolo^-y  (dry) 
Meteorology  (hot) 

Study 

Year 

3 

19^2 

1963 

ilydrolo^’y  (wet) 
Meteorolo/'y  (cold) 

These  three  combinations  of  hydrologic  and  meteorologic  data  were  ana- 
lyzed in  the  initial  simulations.  The  hydrologic  data  for  the  three 
study  years  are  siiown  in  Table  1 . 

58.  Based  on  projected  operations  studies,  mean  monthly  inflows 
and  afterbay  releases  were  provided  by  GPK.  Since  daily  input  is 
required  by  the  mathematical  model , the  mean  monthly  value  was  useci  for 
each  day  of  simulation.  GPK  also  provided  the  daily  operational 
schedules.  Table  2 (generation  and  pumpback  flow  rates  and  durations), 
for  the  three  selected  hydrologic  years  (except  that  tliere  is  no  inunn- 
back  for  Study  Year  2). 

59.  As  discussed  in  the  prelimin;iry  report  (Appendix  A),  a flow 
weighted  average  of  New  Bullard's  Bar  Reservoir  i’ole:isc  flow  temj'era- 
tures  and  South  Yuba  River  temperatures  was  used  to  represent,  t.hi'  in- 
flow temperatures  to  Marysville  Lake  for  all  three  study  ye.-u's. 
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Althouf’h  other  ini'lows  to  the  jjroJect  will  exist,  they  were  not  con- 
sidered in  the  determination  oi‘  iriTlow  temperatures  since  the  New 
Bullard's  Bar  releases  and  South  Yuba  River  flow  constitute  the  majority 
of  tiie  total  project  inflow.  Since  the  preliminary  simulations  showed 
Marysville  Lake  to  be  relatively  inse.nsitive  to  variations  of  input  tem- 
peratures used,  it  was  believed  that  the  use  of  one  set  of  it.flow  temper- 
atures for  all  throe  study  years  was  warranted.  These  flow  te.mperatures 
are  shown  in  Plate  Al. 

60.  Temperature  objectives  at  the  confluence  ol'  tlie  Yuba  and 
Feather  Rivers  were  furnished  by  SPK.  The  objectives  consisted  of  an 
upper  and  lower  limit  on  stream  temp.erature  as  functions  of  time  through 
the  year.  The  objectives  {identified  as  "Downstream  Objectives"  on  the 
various  plates)  are  siiown  on  all  the  plots  of  predicted  temperature  at 
the  release,  afterbay,  and  confluence  ( "Dowistream'' ) . 

Initial  Simulations  and  Model  Calibration 


61 . A preliminary  I'eport  containing  tlie  initial  mathematical  simu- 
lations was  furnished  SPK  in  December  i975-  The  purposes  of  the  initial 
simulations  were  to  (a)  apx.roximate  the  feasibility  of  the  project  to 
meet  downstream  temperature  objectives,  (b)  ietei’mine  port  locations  tha‘ 
allow  tiie  project  to  meet  the  objectives,  (c)  permit  comparison  with  t'ne 
final  numerical  simulations  to  assess  the  benefits  of  the  pipysical  model 
investigation  and  mathematical  model  refinements,  and  (d)  allow  determi- 
nation of  proper  coefficients  by  model  calibration.  That  j'relirainary 
repor't  is  included  as  Appendix  A. 

62 . The  report  of  December  1975  contained  only  simulations  for 
1962  liydrologic  and  meteorologic  conditions.  Subsequently,  fPK  re- 
quested that  initial  simulations  be  performed  foi'  hydroloric  years  193^' 
and  19^*2.  Both  of  these  iiydrologic  years  were  sin;ulated  with  three 
years  of  meteorologic  data,  1962,  1963,  and  1967.  These  simulations 
required  five  selective  witlidrawal  levels  and  a floodgat.e  ti'>  meet 
temperature  objectives.  The  re.sults  wore  furtiished  f.PK  diarini*  Ajril 
1976.  In  May  1976,  CI'K  furni.'.hed  WKf.  with  rovi.-.ed  rout.ings  based  oii  a 
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new  operation  study,  and  the  initial  simulations  were  af^ain  performed. 
Each  of  the  hydrologic  years  193^,  19^2,  and  I962  was  simulated  witii  all 
three  meteorologic  years,  1962,  1963,  and  1967-  These  results  were 
furnished  SPK  during  June  1976. 

63.  As  discussed  previously,  the  matliematical  model  requires  the 
determination  of  coefficients  of  surface  heat  exchange  distribution  and 
internal  mixing.  The  preliminary  report  describes  the  calibration  of 
the  model  for  the  initial  simulations.  The  following  coefficients  were 
determined  from  that  analysis. 


3 = 0.95 
A = 0.3 
= 0.1 
= 0.1 


where 

3 = percentage  of  incoming  shortwave  radiation  absorbed  in  the 
surface  layer 

X = absorption  coefficient 
a,  = mixing  coefficient  at  surface 

X 

= mixing  coefficient  at  bottom 

As  discussed  in  Appendix  A,  the  value  of  3 was  tiecessarily  large  to 
account  for  the  use  of  10-ft  vertical  layers  in  tlie  matiiematical  model. 

. Following  the  initial  simulation.s , tiie  surface  fieat  exciiange 
jortion  of  tiie  matiiematical  model  was  modified  so  tiiat  ttie  3 coef- 
ficient is  applied  to  the  upper  2 ft  of  tiie  lake,  regardless  of  tiie 
vertical  layer  size.  Tiiis  allows  values  for  B and  X tliat  can  be 
more  easily  correlated  to  ligtit  transparency  data.  Additionally,  tiie 
lake  was  divided  into  5-ft  vertical  layers  to  yield  greater  detail  in 
the  model  results.  Tncorporat  i ng  tiiese  ciianges,  tiie  matiiematical  model 
was  recal  iiirat.ed  against  tiie  orig.inal  caliiiration  runs  and  data  for 
1976  on  light,  tr.arispareticy  and  tliormal  profile's  at.  New  Hullard's  Bar 


and  Enclebripht  Reservoirs.  These  data  were  provided  by  GPK.  The  value 
of  the  coefficients  from  this  recalibration  were: 

0.35 
0.12 
0.1 
0.1 

Predicted  thermal  profiles  at  Marysville  Lake  for  1962  hydrologic  and 
meteorologic  data  and  observed  thermal  profiles  at  New  Bullard's  Bar 
Lake  for  1975  and  1976  are  shown  in  Plate  1.  The  predicted  profiles  at 
Marysville  Lake  were  obtained  by  simulating  the  lake  considering  no 
pumpback  flow  and  water  being  released  only  from  the  highest  available 
port.  Comparison  of  the  profiles  in  Plate  1 shows  that  the  predicted 
Marysville  Lake  profiles  have  a similar  shape  and  have  similar  surface 
and  bottom  temperature  ranges  to  the  observed  Now  Bullard's  Bar  Lake 
profiles . 


Final  Simulations 


65.  Mathematical  simulations  were  conducted  for  all  three  study 
years  for  two  pumpback  conditions.  These  conditions  were  pumpback  at 
the  bottom  floodgate  and  pumpback  at  the  level  of  witiuirawal.  These 
simulations  were  first  conducted  with  ttie  me;in  target  temperatures  at 
tiie  confluence  of  tlie  Yuba  and  Feattier  Rivers  as  the  target  releasi' 
temperatures  for  Marysville  Lake.  Tiie  predicted  downstrc:uii  temperatures 
were  not  acceptable,  alUiougii  this  approach  tiad  iieen  used  i!i  the  initial 
simulations  (Appendix  A)  ;ind  yielded  satisfactory  results.  With  the 
revi:!ed  mathematical  model,  the  effect  of  the  afterbay  on  the  downstre;uii 
temfieratures  was  more  real  i stically  (iescribed.  The  effect,  of  tlie  i-e- 
regulation  pool  is.  such  that  downstream  temperature  object.i  ver.  mur.t 
antic  ip.ated . An  optr.ational  change  in  Marysville  Lake  must,  be  mad<> 
with  sufficient  I e.ad  time  to  (iro<luce  .iesired  downstta'ani  t ('mpei'at.utv' 
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at  the  desired  time.  Based  on  the  first  set  of  simulations,  revised 
target  release  temperat\ires  for  Marysville  Lake  were  determined  so  that 
the  resultant  downstream  temperatures  were  within  the  range  of  desired 
downstream  target  temperatures.  The  results  of  all  final  simulations 
are  shown  on  the  following  plates: 

Pumphack  at  the  Bottom 

Predicted  In-Lake  Profiles 
Study  Year  1 - Plate  2 
Study  Year  2 - Plate  3 no  pumphack 
Study  Year  3 - Plate  L 

Predicted  Release,  Afterbay,  and  Downstream  Temperatures 
Study  Year  1 - Plate  5 
Study  Year  2 - Plate  6 no  pumphack 
Study  Year  3 - Plate  7 
Pumphack  at  the  Level  of  V/ithdrawal 

Predicted  In-Lake  Profiles  i 

Study  Year  1 - Plate  8 ^ 

Study  Year  2 - same  as  Plate  3 since  there  is  no  'pumphack  ^ 

Study  Year  3 - Plate  9 \ 

Predicted  Release,  Afterhay,  and  Downstream  Temperatures  i 

Study  Year  1 - Plate  10  1 

Study  Year  2 - same  as  Plate  6 since  there  is  no  pumphack 
Study  Year  3 - Plate  11 

Although  blending  of  flow  from  multiple  vertical  outlets  to  achieve 
temperature  objectives  was  used  in  the  initial  study  and  can  be  quite 
effective,  application  to  a project  such  as  this  could  he  impractical 
from  an  operations  view.  Therefore,  these  simulations  were  conducteii 
assuming  that  releases  from  all  si.x  turbines  occur  from  the  sume  with- 
drawal level  (no  blending). 


PART  V:  DISCUSSION 


66.  The  results  of  the  final  simulations  indicate  that  the  Marys- 
ville Lake  Project  can  meet  temperature  objectives  at  the  Yuba  and 
Feather  Rivers  confluence.  Predicted  downstream  temperatures  were  wit}i- 
in  the  band  of  target  temperatures,  with  only  a few  exceptions,  for  all 
three  study  years  and  for  both  pumpback  conditions  studied.  These  three 
study  years  represent  a broad  range  of  hydrologic  and  meteorologic  con- 
ditions. Because  release  temperatures  can  change  between  Marysville  Lake 
and  the  confluence  of  the  Yuba  and  Feather  Rivers,  different  release 
target  temperatures  had  to  be  established  in  anticipation  of  the  temper- 
ature change.  Release  temperatures  from  Marysville  Lake  can  be  altered 
to  meet  the  downstream  objectives  aL  the  confluence  of  the  Yuba  and 
Feather  Rivero,  as  demonstrated  from  the  fact  that  downstream  tempera- 
tures can  be  moved  about  within  the  objective  temperature  band  witliout 
blending  of  withdrawal  flows.  Tins  again  indicates  the  ability  of  iiie 
project  to  satisfy  objective  temperature  criteria.  At  this  stage  of 

the  project  planning,  it  is  felt  that  to  demonstrate  the  temperatures 
of  the  downstream  releases  would  be  within  an  acceptable  range  should 
be  sufficient.  At  a late’’  stage,  ttie  witlidrawal  structure  design  and 
reservoir  operation  can  be  refined  to  satisfy  more  stringent  downstream 
temperature  objectives,  if  desired. 

67.  Since  botti  pumpback  conditions  stuilied  resulted  in  satisfying 
downstream  temperature  objectives,  the  di f ferf-nce  iti  pumpback  at  tiie 
bottom  or  at  the  level  of  witiidrawal  can  best  be  assessed  by  analysis 
of  the  predicted  in-lake  profiles.  As  could  be  expected,  tlie  effect  of 
pumpback  at  the  bottom  is  predominant  in  the  lower  portion  of  the  lake 
while  pumpback  at  the  level  of  witiidrawal  primarily  al’fected  only  tlie 
upper  levels.  The  pumpback  of  water  at  tiie  bott.om  will  cause  consider- 
able warming  in  the  liypolimnion  and  reduction  of  cold  water  reserves. 

It  should  also  be  noted  that  tiie  results  of  tiie  pliysical  model  indicated 
that  more  mixing  witiiin  tiie  lake  was  induced  by  pumpback  at.  tiie  bottom 
tiian  by  pumpback  at  the  level  of  wi  t iidrawal  . 

68.  A sensitivity  .analysis  w.as  conducte'i  determine  iiow  inucii 
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effect  on  release  temperatures  a significant  error  in  the  entrainment 
coefficient  E would  have  for  pumpback  through  the  floodgate  (at  the 
bottom).  Because  considerable  mixing  occurred  in  the  liypolimnion  for 
this  pumpback  configuration  and  there  was  some  difficulty  in  determin- 
ing E , this  analysis  was  necessary  to  assess  the  reliability  of  tlie 
results.  To  determine  the  impact  of  E on  release  temperatures.  Study 
Year  1 was  simulated  with  all  conditions  (pumpback  at  the  bottom,  port 
selection,  flows,  etc.)  held  constant  except  the  value  of  E which 
was  changed  from  2.5  to  5.0.  The  effects  of  E on  release  temperatures 
from  Marysville  Lake  and  downstream  temperatures  are  presented  in 
Plate  12.  The  dashed-line  circle  represents  a base  condition  of  release 
temperature  computed  by  the  numerical  model  with  E = 2.5  • The  solid 
line  plots  the  deviation  in  temperature  from  the  base  condition  as  a 
result  of  the  perturbation  (E  = 5-0).  The  base  and  perturbed  results 
are  plotted  throughout  the  cycle  of  one  year  (360°  = 1 year)  starting 
on  1 January  positioned  at  the  scaled  axis  and  rotating  counter- 
clockwise. The  results  indicate  that  this  variation  in  E increased 
the  Marysville  Lake  release  temperature  by  a maximum  of  1.5°C  and  the  ) 

downstream  temperature  by  a maximum  of  0.75°C  (Plate  12).  The  temper- 
ature increase  was  realized  only  during  the  last  quarter  of  the  year 
when  releases  were  made  from  the  floodgate.  The  small  variation  in 
downstream  temperatures  indicates  that  the  downstream  temperature 
routings  are  not  very  sensitive  to  E and  that  confidence  in  the  pre- 
diction can  be  maintained  even  with  a large  error  in  the  E value. 

69.  The  piiysical  model  provided  valuable  understanding  of  Uie 
mixing  characteristics  within  Marysville  Lake  and  afterhay.  Results 
from  the  physical  model  tests  were  used  to  modify  the  mathematical  model. 

These  modifications  described  the  effects  of  the  generation  and  pumpback 
operations  on  the  density  structure  within  Mai’ysville  Lake.  Result;-, 
were  also  used  to  describe  the  proportions  of  previous  generation  water 
and  afterbay  water  that  composed  the  pumpback  flow.  Tlii;-.  al  lowiai  proper 
simulation  of  tiie  pumpback  temperature  and  mixing  of  re;-.idu.-il  w.-iter  in 
the  afterbay . 

70.  The  results  of  the  final  simulations  for  I'h-.’  (I’l.ates  2 .-md  ''■) 
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with  pumpback  at  the  bottom  can  be  compared  with  the  initial  simulation 
results  shown  in  Appendix  A.  The  predicted  thermal  profiles  within  the 
lake  are  similar  even  though  the  entrainment  by  the  pumpback  Jet  was 
2.5  times  greater  in  the  final  simulations  than  it  was  in  the  initial 
simulations.  This  is  because  in  both  cases  the  pumpback  current  was 
placed  into  the  region  of  the  lake  below  the  thermocline  where  tempera- 
ture differences  are  small  and  changes  in  the  thermal  structure  would 
not  be  great.  Although  the  final  and  initial  predicted  downstream 
temperatures  are  similar,  the  target  release  temperatures  for  Marysville 
Lake  to  achieve  those  downstream  temperatures  were  quite  different. 

Also  the  release  and  afterbay  temperatures  are  different  in  the  initial 
and  final  simulations.  In  the  final  simulations,  the  effect  of  the 
afterbay  is  properly  accounted  for,  and  the  results  show  that  Marysville 
Lake  must  be  operated  to  complement  or  compensate  for  the  afterbay 
effect  (paragraph  65)  in  order  to  achieve  the  desired  downstream  temper- 
atures. Certainly,  the  final  mathematical  model  more  realistically 
describes  the  physical  situation  than  the  initial  mathematical  model  did 
and  should  serve  as  an  excellent  tool  for  design  of  the  withdrawal 
structure . 

71.  All  port  locations  were  used  by  the  mathematical  model  during 
the  final  simulations.  The  port  operation  schemes  used  in  the  simu- 
lations are  described  at  the  top  of  the  predicted  release  temperature 
plots  of  Plates  5?  6,  7,  10,  and  11.  A summary  of  the  ports  used  for 
each  study  year  and  pumpback  condition  follows: 


Study  Year 

1 1962  Hydrology 

1962  Meteorology 

2 193^  Hydrology 
1967  Meteorology 

3 19^*2  Hydrology 

1963  Meteorology 


Location  of  Pumpback 

Withdrawal  level 
Bottom 

No  pumpback 


Ports  Used 

2,  3,  6 
2,  3,  6 

5,  6 


Withdrawal  level  1,  2,  h* , 5*,  6 

Bottom  1,  2,  6 


* Ports  I4  and  5 were  used  for  very  short  periods  and  were  not 
necessary  for  this  case. 
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The  port  numbers  represent  the  following  port  locations: 


Port 
1 
2 
3 
h 

5 

6 (floodgate) 

From  these  results  it  appears  that  port  5 could  possibly  be  eliminated. 

At  least  four  ports  and  a floodgate  were  necessary  because  of  the  large 
fluctuation  in  pool  level  (fluctuates  from  el  h29  to  56o  over  the  three 
years  studied)  that  occurs  over  the  three  study  years.  Witii  less  fluctu- 
ation in  pool  level,  the  temperature  objectives  could  be  satisfied  with 
fewer  j)orts.  Throughout  this  study,  port  6 has  been  referred  to  as  a 
floodgate.  Because  the  flow  capacity  of  all  ports  was  sufficiently 
large  to  handle  all  outflows  simulated,  it  was  not  necessary  to  operate 
port  6 as  a floodgate  to  release  large  flows.  Therefore,  port  6 can  be 
considered  a low  level  withdrawal  port  rather  than  a floodgate. 

72.  It  should  be  emphasized  tliat  the  port  locations  used  in  this 
study  are  not  tiie  only  outlet  locations  wiiicii  could  be  used  to  satisfy 
temperature  objectives.  For  example,  with  a pool  level  that  fluctuates 
substantially,  it  could  be  difi'icult  at  times  to  meet  warm  water  ob- 
.Jectives  with  fixed  ports.  Movable  ports  or  weirs  may  be  more  effi'ctive 
in  aciiieving  temperature  objectives.  The  matiiemat ical  model  in  con- 
junction witii  tile  ptiysical  model  can  be  used  in  tlie  design  pii.ase  of  the 
project  to  evaluate  ports,  weirs,  and  combinations  of  ports  and  weirs 
that  would  satisfy  the  temfierature  objectives.  ^election  of  the  final 
liesign  could  tiien  l>c  made  considering  economics,  tiy.iraul ics , and  opei'- 
ational  factors. 
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Table  1 

Hydrolopiic  Input  Data,  Hydrolo^uc  Year^; 
IQ62,  193^4,  and 


Year 

1962 


1934 


19^*2 


Month 

Inflow, 
cfs  per  day 

Afterbay 
Release , 
cfs  per  day 

Irrigat ion 
Withidrawal  at 
Daguerre  Point, 
cfs  per  Jay 

Jan 

l!4l49 

II43I4 

0 

Feb 

338i* 

1798 

0 

Mar 

22kJ 

IOI45 

50 

Apr 

2237 

2172 

310 

May 

2136 

2182 

800 

Jun 

236U 

2I42I4 

1000 

Jul 

2288 

2273 

1130 

Aug 

17TA 

I8I43 

1030 

Sep 

1283 

I3I43 

800 

Oct 

i*89i4 

1*975 

360 

Nov 

1010 

2020 

0 

Dec 

2061 

1985 

0 

Jan 

1773 

1758 

0 

Feb 

2081 

2076 

0 

Mar 

192i* 

I92I4 

50 

Apr 

lit  3^ 

1150 

31*0 

May 

1323 

IU85 

800 

Jun 

1308 

131*8 

8U0 

Jul 

lll4i 

II4I4I4 

Q30 

Aug 

1051 

1371* 

860 

Sep 

828 

1187 

690 

Oct 

753 

1020 

31*0 

Nov 

707 

707 

0 

Dec 

631 

682 

0 

Jan 

5722 

5551 

0 

Feb 

6758 

639I* 

0 

Mar 

2556 

1111 

50 

Apr 

^*793 

292I4 

310 

May 

14672 

I4lt0l4 

800 

Jun 

3788 

3788 

1000 

.Jul 

■091 

3056 

1130 

Aug 

1717 

1793 

1130 

Sep 

1258 

2929 

800 

Oct 

995 

2OI4O 

liOO 

Nov 

71*2 

1753 

0 

Dec 

2177 

2096 

0 

I 

I 
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Table  2 


Year 

1962 


I93U 


19^2 


Operation  Schedule,  Hydrologic 

Years 

1962,  1931+ , and  I9U2 

Generation 

Pumpback 

Generation 

Duration, 

Pumpback 

Duration, 

Month 

Rate,  cfs 

hr  per  day 

Rate,  cfs 

hr  per  day 

Jan 

39,800 

1.68 

2960 

9.00 

Feb 

38,900 

1.88 

2560 

9.00 

Mar 

36,900 

1.73 

3830 

9.00 

Apr 

36,i+00 

1.80 

1000 

2.1+2 

May 

36,500 

I.7I+ 

1000 

2.1+2 

Jun 

36,900 

2.33 

2020 

9.00 

Jul 

37,700 

3.71 

8200 

9.00 

Aug 

37,500 

3.71+ 

9860 

9.00 

Sep 

37,^400 

3.28 

9280 

9.00 

Oct 

38,300 

3.63 

0 

0 

Nov 

37,500 

1.78 

lll+O 

9.00 

Dec 

38,100 

1-71 

1000 

7.62 

Jan 

U3,200 

I.1I+ 

No  pumpback 

for  this 

Feb 

i*3,000 

1.35 

study  ; 

year 

Mar 

i*3,200 

1.25 

Apr 

1+3,500 

0.93 

May 

1+3,500 

0.95 

Jun 

1+3,000 

0.88 

Jul 

1+2,600 

0.95 

Aug 

1+2,100 

0.91 

Sep 

1+1,300 

0.80 

Oct 

1+1,000 

0.70 

Nov 

1+0,800 

0.1+8 

Dec 

1+0,800 

0.1+7 

Jan 

39,800 

3.90 

0 

0 

Feb 

39,900 

1+.1+8 

0 

0 

Mar 

36,500 

1.73 

3580 

9.00 

Apr 

35,1+00 

2.1+0 

0 

0 

May 

3l+,900 

3.53 

0 

0 

Jun 

3l+,500 

3.07 

0 

0 

.lul 

3l+,800 

3.83 

I+9I+O 

9.00 

Aug 

3l+,700 

3.82 

9170 

9.00 

Sep 

35,700 

3.33 

3750 

9.00 

Oct 

36,500 

2.26 

281+0 

9.00 

Nov 

37,500 

1.78 

1980 

Q.OO 

Dec 

38,100 

1.71 

1000 

1+.1+8 

7/17/62- 
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APPENDIX  A:  MARYSVILLE  LAKE  TEMPERATURE  STUDY 

INITIAL  MATHEMATICAL  MODEL  INVESTIGATIONS 

PART  I : INTRODUCTION 

1.  The  Marysvi.lle  Lake  Temperature  Study  is  currently  being  con- 
ducted at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  for 
the  U.  S.  Army  Engineer  District  Sacramento  (SPK).  Both  a physical 
model  and  a numerical  simulation  model  will  be  used  during  the  study. 

The  study  will  consist  of  an  initial  mathematical  model  investigation, 
prior  to  the  construction  of  the  physical  model,  and  then  a detailed 
investigation  using  both  physical  and  mathematical  models.  This  pro- 
cedure was  used  (a)  to  provide  an  early  insight  into  the  ability  of 
Marysville  Lake  to  meet  downstream  temperature  objectives,  (b)  to  permit 
comparison  of  results  of  the  numerical  simulations  with  and  without  the 
iiydrodynamics  as  defined  by  the  physical  model,  thereby  providing  some 
assessment  of  the  benefit  of  a physical  model  investigation,  and  (c)  to 
allow  determination  of  proper  values  of  coefficients  to  be  used  in  sub- 
sequent simulations  in  order  to  provide  the  most  reliable  results  possi- 
ble from  tiie  simulations.  The  purpose  of  this  preliminary  report  is  to 
present  the  results  of  the  initial  mathematical  model  investigations. 

2.  The  proposed  Marysville  Lake,  witli  the  dam  at  the  Parks  Bar 
Site  on  the  Yuba  River,  will  be  about  360  ft  deep  witli  a volume  of 
916,000  acre- ft  imd  an  area  of  6,6hO  acres  at  gross  pool  el  560.  A 
selective-withdrawal  facility  will  provide  flows  to  the  powerhouse.  The 
project  provides  for  pumped-storage  power  operations,  using  an  afterbay 
immediately  downstream  of  the  dam. 

3.  The  afterbay  of  Marysville  Lake  will  be  about  73  ft  deep  with  a 
volume  of  )tU,000  acre-ft  and  an  area  of  9^*0  acres  at  gross  pool  el  233* 

A gated  overflow  spillway  with  sluices  will  control  flow  from  the  after- 
bay. 

It.  The  Yuba  River  joins  tiie  I’eather  River  approximately  lU  miles 
b(;low  the  Marysville  D;un.  Tt  is  at  this  confluence  that  temperature- 
objectives  arc  estalil  irlicd . 
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PAKT  II;  MODEL  DESCRIPTION 


5.  The  dovrtu:tre!un  relea:;e  characteristics  and  the  internal  struc- 
ture of  temperature  for  Mar’yrville  Lake  were  predicted  usinR  a numerical 
simulation  model.  Ilie  model  (WEDTPiX)  used  in  conjunction  with  this 

invest ir.a-tion  wa;;  developed  at  WES  based  upon  the  results  of  Clay  and 

I*  2 h 5 

Fruh,  Edinger  and  Geyer,  Dake  and  Harleman,  and  Bohan  and  Grace. 

6.  The  WESTEX  model  provides  a procedure  for  examining  the  balance 
of  thermal  energy  imposed  on  an  impoiindment . These  energy  balance  and 
lake  hydrodynamic  phenomena  are  used  to  map  vertical  profiles  of  temper- 
ature in  the  time  domain.  The  model  includes  computational  methods  for 
simiHating  heat  transfer  at  the  air-water  interface,  heat  advection  due 
to  inflow  and  outflow,  and  the  internal  dispersion  of  thermal  energy. 

The  model  is  conceptually  based  on  the  division  of  the  impoundment  into 
discrete  horizontal  layers.  Fundamental  assumptions  include  the  follow- 
ing: 

Isotherms  are  parallel  to  the  water  surface  both  laterally 
and  longitudinally. 

b.  Tiie  water  in  each  discrete  layer  is  isotropic  and  physi- 
cally homogeneous. 

£.  Internal  advection  (between  layers  within  the  lake)  and 
heat  transfer  occur  only  in  the  vertical  direction. 

External  advection  (inflow  to  and  outflow  from  the  lake) 
occurs  as  a uniform  horizontal  distribution  within  each 
layer. 

e.  Internal  dispersion  of  thermal  energy  is  accomplished  by 
a diffusion  mechanism  that  combines  the  effects  of  molecu- 
lar diffusion,  turbulent  diffusion,  and  thermal  convection. 

7.  Tlie  surface  heat  exchange,  internal  mixing,  inflow,  and  out.- 
flow  processes  are  simulated  separately,  and  their  effects  are  intro- 
duced sequentially  at  daily  intervals. 

fi.  The  WEGTEX  model  employs  an  a]iproach  to  the  evaluation  of 

2 

surface  heat  transfer  develo[jed  by  Edinger  and  Geyer.  Tliis  method 
formulates  equilibrium  temperatures  and  coefficients  of  surface  heat 


* Gee  References,  pa/;e  All. 


A? 


F 


1 


exchanf^e.  Equilibrium  temperat\u'e  is  defined  as  that  temperature  at 
which  the  net  rate  of  heat  exchange  between  the  water  surface  and  the 
atmosphere  is  zero.  The  coefficient  of  surface  heat  exchange  is  the 
rate  at  which  the  heat  transfer  process  occurs.  The  equation  describ- 
ing this  relationship  is: 


11  = K(E  - 9) 


(1) 


where 

2 

H = net  rate  of  surface  heat  transfer,  Btu/ft  /day 
K = coefficient  of  surface  heat  exchange,  Btu/ft^/ day /°F 
E = equilibrium  temperature , °F 
9 = surface  temperature,  °F 

The  computation  of  equilibrium  temperature  and  heat  exchange  coefficient 
is  based  solely  on  meteorological  data  as  outlined  by  Edinger,  Duttweiler, 

3 

and  Geyer. 

9.  The  net  heat  exchange  at  the  surface  is  composed  of  seven  heat 
exchange  processes : 

Shortwave  solar  radiation, 
b.  Reflected  shortwave  radiation. 

£.  Long-wave  atmospheric  radiation. 

Reflected  long-wave  radiation. 

£.  Heat  transfer  due  to  conduction. 

_f.  Back  radiation  from  the  water  surface. 

£.  Heat  loss  due  to  evaporation. 

For  every  day  of  meteorological  data,  the  seven  heat  exchange  terms  can 
be  evaluated  and  the  not  heat  exchange  rate  expressed  in  term::  of  an 
equilibriiun  temperature  and  an  exchruige  coefficient. 

10.  All  of  the  surface  heat  exchange  pi-ocesses,  with  the  exception 
of  shortwave  radiation,  affect  only  the  top  few  feet  of  the  lake.  Sdiort.- 
wave  radiation  penetrates  the  water  siurface  and  increases  the  temper- 
ature at  greater  depths.  Rased  on  laboratory  investigations,  Dakc  and 
Harleman  have  su/igested  an  exponential  decay  with  depth  for  describing 
the  heat  flux  due  to  shortwave  penetration. 
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11.  The  surface  }ieat  exchange  concepts  are  implemented  in  the 
'ifESTEX  model  by  tlie  exiionential  penetration  of  a percentage  of  the  in- 
coming shortwave  radiation  and  the  placement  of  the  effect  of  all  other 
sources  of  surface  heat  exchange  into  the  surface  layer.  This  can  be 
expressed  mathematically  by  the  following  two  equations: 


Hg  = K(E  - Q)  - (1  - 6)S 


(2) 


H.  = (1  - 6)  Se  ^^i 


(3) 


where 

6 

S 

H. 

1 

e 

X 

z . 
1 


2 

heat  transfer  rate  into  or  out  of  surface  layer,  Btu/ft  /day 

shortwave  radiation  absorbed  in  the  surface  layer,  percent 

2 

rate  of  total  incoming  shortwave  radiation,  Btu/ft  /day 

2 

rate  of  heat  absorption  in  layer  (i),  Btu/ft  /day 
natural  logarithmic  base  (2.7183) 
absorption  coefficient,  ft  ^ 
depth  below  surface,  ft 
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12.  Tlie  process  of  inflow  into  a lake  is  simulated  in  WESTEX  by 
the  placement  of  Inflow  quantity  and  quality  at  that  layer  in  whicli  the 
density  of  the  lake  corresponds  most  nearly  to  the  density  of  the  in- 
flow. Research  efforts  and  physical  model  studies  at  VIES  indicate  the 
existence  of  entrainment-induced  density  currents  which  flow  upstretun 
along  the  surface  into  the  turbulent  mixing  zone  caused  by  the  inflow. 
Entrainment  is  implemented  in  the  model  by  aufqnenting  the  inflow  quan- 
tity with  a volume  from  the  surface  layer  (as  discussed  in  paragraph  13). 
Characteristics  of  inflow  and  the  entrained  flow  are  averaged,  and  mixed 
values  of  density,  temperature,  and  other  water-quality  parameters  are 
determined.  The  mixed  density  is  used  to  determine  placement  of  the 
total  quantity  .'and  mixed  quality.  In  the  model  the  process  of  inflow 
displaces  upward  a volume  equal  to  the  total  inflow  plus  entrainment 
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quajitity.  This  upward  displacement  is  reflected  in  the  model  by  an 
increase  in  the  water  surface.  A corresponding  decrease  in  water  surface 
occurs  as  a result  of  the  outflow  process. 

13.  The  volume  of  the  entrained  current  is  generally  expressed  as 
a percentage  of  the  inflow  quantity.  Prior  flume  studies  indicated  that 
this  percentage  ranges  from  25  to  nearly  200.  The  percentage  is  thought 
to  be  a function  of  slope,  width,  flow  quantity,  density  of  inflow,  and 
density  within  the  lake,  but  analytical  relationships  have  not  yet  been 
determined.  Thus,  a physical  model  is  of  significant  benefit  in  the 
evaluation  of  entrainment.  For  these  initial  simulations,  it  was  as- 
sumed that  there  is  no  entrainment  current  generated  by  the  inflow. 

ll.  The  inti,  rnal  mixing  process  is  simulated  by  a transfer  of 
heat  and  other  water-quality  constituents  between  adjacent  layers.  The 
magnitude  of  the  heat  transfer  between  two  layers  is  a percentage  of 
the  total  heat  transfer  required  to  completely  mix  the  two  layers.  The 
percentage  is  a mixing  coefficient  that  is  definea  for  every  layer. 

Data  input  includes  values  of  the  mixing  coefficient  at  the  top  and  at 
the  bottom  of  the  lake.  An  exponential  fit  between  the  two  extreme 
values  is  used  to  determine  the  appropriate  coefficient  at  each  layer. 
Because  of  the  lack  of  knowledge  concerning  diffusion  processes  and  the 
lack  of  Justification  for  a more  complex  procedure,  this  simple  dif- 
fusion Einalogy  is  applied  to  mixing  between  layers  for  temperature  and 
other  water-quality  parameters. 

15-  The  outflow  component  of  the  model  incorporates  the  selective- 
withdrawal  techniques  developed  at  WES.^  Transcendental  equations  de- 
fining the  zero  velocity  limits  of  the  withdrawal  zone  are  solved  with 
a half-interval  search  method.  With  knowledge  of  the  withdrawal  limits, 
the  velocity  profile  due  to  outflow  can  be  determined.  The  flow  from 
each  layer  is  then  the  product  of  the  velocity  in  the  layer,  the  width 
of  the  layer,  and  the  thickness  of  the  layer.  A flow-weighted  average 
is  applied  to  water-quality  profiles  to  determine  the  value  of  the  re- 
lease content  of  each  par^imetor  for  each  time  step. 

16.  Tlie  lake  regulation  algorithms  have  been  developed  to  real- 
istically simulate  the  field  operation  of  a selective-withdrawal  system. 
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The  selective-withdrawal  system  is  assijmed  to  he  configured  with  an 
arbitrary  number  of  selective-withdrawal  intakes  located  in  each  of  two 
wet  wells  with  a separate  floodgate.  Maximum  and  minimum  flows  from 
each  intake  and  from  the  floodgate  must  be  specified.  Also,  the  maximum 
flow  for  the  selective-withdrawal  system  is  specified.  The  algorithms 
attempt  to  numerically  withdraw  water  at  or  near  the  objective  temper- 
ature. Withdrawal  will  be  from  either  one  intake  level,  two  adjacent 
intake  levels,  and/or  the  flood-control  intake,  depending  upon  the  ob- 
jective temperature,  the  temperature  profile,  the  int-ike  capabilities, 
and  the  amount  of  flow  to  be  released. 

IT.  humpback  currents  are  represented  in  the  model  with  the  three 
processes  of  entrainment,  mixing,  and  placement.  Entrainment  is  ex- 
pressed as  a percentage  of  pumphack  volume.  The  entrained  water  at  its 
temperature  is  numerically  mixed  with  the  pumpback  volume  at  its  temper- 
ature. This  total  volume  of  water  at  the  mixed  temperature  is  placed 
into  the  lake  at  the  appropriate  density  level.  The  vertical  location 
from  which  entrainment  takes  place  and  the  percentage  of  entrainment  are 
required  as  data  input  to  the  mode] . A physical  model  is  of  great  bene- 
fit for  assessment  of  these  parameters. 
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PART  III:  DEVELOPMENT  OF  INPUT  DATA 


l8.  At  the  request  of  SPK  an  average  hydrologic  year  (1962)  was 
used  for  the  initial  simulations.  For  calibration  purposes  simulations 
were  conducted  with  meteorological  data  from  197^  and  1975  and  hydro- 
logical  data  from  1962. 

19-  Meteorological  data  from  Beale  Air  Force  Base  were  used  for 
this  study.  The  required  data  consisted  of  dry  bulb  temperature,  dew 
point  temperature,  wind  speed,  and  cloud  cover.  For  1962,  these  data 
were  obtained  from  the  National  Climatic  Center,  Asheville,  North 
Carolina.  For  197^  Mid  1975,  these  data  were  obtained  from  SPK. 

20.  Based  on  a projected  operation  study,  mean  monthly  Inflows 
and  outflows  were  provided  by  SPK.  Since  daily  input  is  required  by 
the  WESTEX  model,  the  mean  monthly  value  was  used  for  each  day  of  the 
appropriate  month. 

21.  Although  New  Bullard's  Bar  Reservoir  is  immediately  upstream 
of  Marysville  Lake  and  will  contribute  the  majority  of  inflow  to  the 
project.  New  Bullard's  Bar  Reservoir  was  not  numerically  simulated. 

This  decision  was  based  on  the  lack  of  adequate  input  data  for  a numeri- 
cal simulation  as  well  as  uncertainties  about  the  selective  withdrawal 
chiaracteristics  of  the  New  Bullard's  Bar  Intake  structure.  Observed 
stream  temperatui’e  data  were  available  below  Colgate  poweriiou.se,  the 
location  at  wiiicli  t!ie  New  Bullard's  Bar  releases  are  dischai'ged  to  the 
Yuba  River.  By  combining  these  data  with  observed  stream  temperature 
data  on  the  upstre;un  Couth  Y ;ba  River,  a range  of  input  stre.'im  temper- 
ature data  was  established.  Tiiree  sets  of  input  stream  temperatures 
were  tested  for  the  initial  simulations.  The  results  of  the  simulations 
were  relatively  insensitive  to  the  variations  of  input  temperatures 
used;  tiierefore,  one  set  of  input  stream  temperatures  was  used  for  all 
subsequent  simulations.  Plate  A1  shows  the  inflow  temperatures  tiiat 
were  used.  These  temperatures  were  computed  as  a flow-weighted  average 
of  New  Bullard's  Bar  Reservoir  release  flows  aril  South  Yuba  River  flows. 

22.  Temperature  ob.lectives  at  tlic  confluence  of  the  Yuba  and 
Featiier  Rivers  wore  furn i I'.lie.i  to  WE2.  by  SPK.  The  objectives  consisted 
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PART  IV:  MODEL  CALIBRATION 


23.  Aa  has  been  discussed  previously,  the  WESTEX  model  requires 
the  determination  of  coefficients  of  surface  heat  exchange  distribution 
and  internal  mixing.  For  Marysville  Lake,  these  coefficients  were  de- 
termined by  conducting  simulations  with  197^  and  1975  meteorologic  data. 
Coefficients  were  adjusted  and  the  simulation  was  repeated  until  the 
predicted  temperature  profiles  corresponded  in  shape  and  range  to  those 
observed  during  1975  in  the  nearby  New  Bullard's  Bar  and  Englebright 
Reservoirs.  The  operation  of  the  outlet  works  used  in  these  simulations 
was  representative  of  the  New  Bullard's  Bar  Reservoir.  The  following 
coefficients  were  determined  from  this  analysis; 

B = 0.95 
X = 0.3 
= 0.1 
«2  = 0-1 


where 

B = percentage  of  incoming  shortwave  radiation  absorlod  in  the 
surface  layer 

X = absorption  coefficient 

= mixing  coefficient  at  surface 
= mixing  coefficient  at  bottom 

2h . It  should  be  noted  that  tiie  value  selected  for  B (0.95)  is 
large  relative  to  values  used  by  WEE  engineers  in  previous  thennal 
studios.  A large  value  of  6 might  traditionally  be  used  to  reprer.c -n' 
an  extremely  turbid  impoundment,  that  is,  one  in  which  light  penetra- 
tions were  sev('roly  restricted.  In  ttio  past,  6 tins  been  selecteu  for 
a 2-ft-laycr  ttiickness  and,  altiiougti  it  is  assumed  tliat  Marys. vi  1 1 c L.aki' 
will  be  relatively  clear,  tliis  v.alue  of  B was  netaied  to  accoimt  fi>r 
the  use  of  10- ft  vortical  layers  in  the  math«'mat,  i c.al  simulations. 
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PART  V;  MODIFICATIONS 
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25-  Several  modifications  were  made  to  the  WESTEX  model  to  acco'orit 
for  the  effects  of  pumped  storage  operations  at  Marysville  Lake  and  to 
compute  the  stream  temperature  at  the  confluence  of  the  Yuba  and  Feather 
Rivers . 

26.  The  MESTEX  model  is  incapable  of  similLating  exactly  the  unique 
configuration  of  the  intake  structure  proposed  for  Marysville  Lake.  The 
model  requires  that  flows  be  released  through  one  or  two  wet-wells  con- 
taining selective-withdrawal  intakes  and  a single  floodgate.  The  pro- 
posed structure  consists  of  six  wet  wells,  each  with  a floodgate  and 
10-ft  temperature  shutters  from  just  above  the  floodgate  to  tlie  top  of 
the  tower.  For  the  initial  simulations,  one  floodgate  level  and  three 
levels  of  selective-withdrawal  intakes  were  assumed.  The  center-line 
elevation  of  the  floodgate  level  was  250  ft  iuid  the  center-line  ele- 
vations of  the  select! ve-withdrawal  intakes  were  375,  ^*'*0,  and  500  ft. 
Blending  of  flow  was  allowed  by  combining  tiie  six  wet  wells  itito  two. 
This  pi’ovided  for  blending  between  adjacent  j'ort  levels  and/or  t!ie 
floodgate  wiiicii  greatly  reduced  trie  numerous  possibilities  by  w:.icn 
blending  could  bo  actiieved  witti  six  wet  w.jlls.  All  punphack  flow.-  wcr<; 
assiuiied  to  pass  throug.h  the  floodgate. 

27.  To  enable  tiie  WEOTEX  model  to  account  for  pumped  storage 
operations,  it  was  necessary  to  account  for  the  effect  of  the  afterb:iy 
on  tile  temperature  of  tlie  generated  release.  It  was  assumed  tliat  tlie 
afteriiay  will  remain  fully  mixed.  Tlie  lieat  excliange  at  the  air-water 
interface  in  tlie  afterbay  wan  taken  into  account  in  a manner  similar  to 
that  for  Marysville  Lake.  Tlie  combined  effect  of  this  heat  exchange  at 
the  air-water  Interface  :uid  the  temperature  of  the  generation  flows 
yieldeil  the  ii.ixed  temperature  of  the  afterbay.  This  tempora+.ure  war. 
used  an  the  temperature  of  the  piunpback  flow  and  of  the  afterbay  release. 

28.  An  observation  of  the  Richard  B.  Russell  Lake  Wat.er  Quality 
Ftudy  * w.as  that  a pumpback  jot  entering  the  reservoir  entrained  an 
lunount  of  reservoir  watei’  resulting  in  'ui  upstream  current  greater  in 
magnitude  than  the  -unount  of  flow  bc'ing  (nimped.  For  M.arysvillc'  wilhout. 
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the  benefit  of  the  physical  model,  this  entrained  flow  was  assumed  to 
come  from  the  area  of  the  reservoir  immediately  above  the  floodgate 
opening,  'fwo  values  of  entrainment  were  assumed,  0 and  100  percent. 

29.  A daily  power  operation  schedule  has  not  been  established  for 

Marysville  Lake.  For  this  study,  GPK  provided  WES  with  maximum  and 

minimum  generation  flow  rates.  Assuming  a i^-hr  generation  period,  the 
generation  flow  rates  analyzed  were  1 ,500,  2^1,000,  and  50,000  cfs.  As 
a base  condition,  one  simulation  1 as  made  assuming  no  generation  or 
pumpback  flow. 

30.  Since  the  objective  temperatures  for  this  project  are  at  the 

confluence  of  the  Yuba  and  Feather  Rivers,  the  afterbay  release  temper- 

atures must  be  routed  approximately  10  miles  downstreain.  This  routing 

is  accomplished  by  assuming  a fully  mixed  flow  and  using  the  relation- 

2 

ship  of  Edinger  and  Geyer. 

T = (T„  - E)  e“  + E 
X 0 

a = -Kx/yC^yU 


where 

T 

X 

T 

0 

E 

K 

X 

Y 

C 

P 

y 

u 


downstre:un  water  temperature,  °F 
initial  water  temperature,  °F 
equilibrium  temperature,  °F 

coefficient  of  surface  lieat  exciiange,  Rtu/ft^/day/'^F 

downst.re.'uii  distance,  ft 

specific  weigiit  ol’  water,  62.^1  Ib/ft^ 

specific  iieat.  of  water,  1.0  Btu/lb-°F 

mean  stream  deptii,  ft 

mean  .-.trerun  vi-locity,  ft/day 


Ttic  values  of  »■  qui  L i in- i un'  ti-mperatur<>  .'ind  i;oefficient  of  s.urface  Innat 
exciiangi'  were  t.hose  u.'.ed  for  Marys.vllle  Lake  foi'  tlie  approprial-e  day. 


All 


Additionally,  for  the  initial  simulations,  no  attempt  was  made  to 
incrementally  route  flow  down  the  Yuba,  but  instead  a constant  value  of 
mean  stream  depth  (6  ft)  and  mean  streeun  velocity  (6  fps)  were  used  to 
accomplish  the  routing  in  one  step. 
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PART  VT:  DISCUSSION 


31.  The  results  of  the  simulations  in  the  form  of  time-histories 
of  computed  temperature  of  the  release,  in  the  afterbay,  and  downstream 
for  1962  and  computed  thermal  profiles  in  Marysville  Lake  for  various 
dates  durine  the  period  April-September  1962  are  presented  in  Plates  A2- 
A7.  Plates  A2  and  A3  present  the  base  condition  results  with  no  power 
flow.  Plates  Ah  and  A5  present  the  results  for  generation  and  pumpback 
flows  with  the  assumption  of  100  percent  entrainment  of  flow  by  the 
pumpback  Jet.  Plates  A6  and  AT  present  the  results  for  generation  and 
pumpback  flows  assuming  no  entrainment  of  flow  by  the  pumpback  jet. 

32.  The  results  indicate  the  ability  of  the  project  to  meet  temper- 
ature objectives.  Inspection  of  the  computed  Marysville  Lake  thermal 
profiles  shows  the  lake  to  be  less  thermally  stratified  for  the  gener- 
ation and  pumpback  flows  with  100  percent  pumpback  entrainment  than 

flow  with  no  entrainment.  It  is  believed  that  the  ability  of  the  pro- 
ject to  meet  temperature  objectives  is  primarily  related  to  the  hydro- 
dynrunics  of  the  lake.  Caution  should  therefore  be  exercised  in  drawing 
any  conclusions  from  these  results,  since  simplifying  assumptions  about 
the  hydrodynamics  were  made  for  this  iuialysis.  The  results  of  the 
physical  and  advanced  mathematical  modeling  efforts  on  the  Marysville 
Lake  Project  will  determine  the  validity  of  these  assumptions.  Only 
after  the  hydrodynamics  of  Marysville  L;ike  have  been  properly  defined 
can  any  modifications  to  the  project  be  analyzed. 
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APPENDIX  B;  SELECTIVE  WITHDRAWAL  ANALYSIS 


3* 

I.  The  WES  generalized  selective  withdrawal  teciinique  uses  the 
following  equation  to  compute  the  upper  and  lower  limits  of  withdrawal 
for  an  orifice: 


Z 


(Bl) 


where 

Z = vertical  distance  from  the  elevation  of  the  orifice  center 
line  (^)  to  the  upper  or  lower  limit  of  the  zone  of  with- 
drawal , ft 

V^  = average  velocity  through  the  orifice,  fps 

A^  = area  of  the  orifice  opening,  sq  ft 

Ap  = density  difference  of  fluid  between  the  elevations  of  the 

orifice  ^ and  the  upper  or  lower  limit  of  the  zone  of  with- 
drawal , g/cc 

= fluid  density  at  the  elevation  of  the  orifice  g/cc 

2 

g = acceleration  due  to  gravity,  ft/sec 
The  product  of  and  A^  gives  the  outflow  rate  (in  cfs)  through 

the  orifice.  Tiierefore,  knowing  the  port  discharge,  Z can  be  deiern.itied 
from  Equat ion  Bl . 

3 

2.  The  selective  withdrawal  studies  were  conducted,  for  single  and 
multiple  vertical  orifices,  so  tiiere  is  no  problem  in  applying  tlie  tecti- 
nique  to  a single  horizontal  port.  However,  wlien  applying  ttie  teciinique 
to  simultfuieous  discharge  througli  multiple  horizontal  ports  (all  at  tiie 
s;ime  level),  the  question  has  been  raised  wtietlier  t.o  use  in  computing 
the  wit.hilrawal  limits  the  (iisciiarge  •uid  port  area  for  a single  port  or 
tile  sum  of  the  discliarg.es  and  cumulal.ive  port  area  for  all  tlie  ports 
operating  at  tiiat  level.  Because  the  preliminary  design  for  t.he 


* gee  He  ferencer  , page  . 


Bl 


Marysville  project  includes  six  ports  at  each  level,  this  question  had 
to  be  answered. 

3.  The  logical  application  of  Equation  B1  was  to  use  the  discharge 
and  port  area  for  a single  port  since  the  equation  war.  developed  for  a 
single  port.  Assioming  that  equal  simultaneous  releases  would  occur 
through  each  port,  the  withdrawal  characteristics  and  release  temper- 
atures for  a single  port  would  apply  to  all  the  ports  at  that  level. 

The  validity  of  this  approach  was  tested  before  incorporation  into  tlie 
Marysville  simulation  model. 

A series  of  tests  were  conducted  in  an  existing  3-ft-wide  by 
2-ft-deep  by  25-ft-long  flume  constructed  of  transparent  plastic  to 
evaluate  the  withdrawal  characteristics  of  multiple  ports  releasing 
simultaneously  from  the  same  level  of  a density  stratified  pool.  Be- 
cause of  the  dimensional  limitation  of  the  flume  widtli,  it  was  only 
possible  to  model  four  of  the  six  ports.  This  was  considered  to  be 
sufficient  for  determining  the  effect  of  multiple  horizontal  ports. 

Each  port  was  0.1  ft  wide  by  0.02  ft  high  or  had  a width  five  times 
larger  than  the  height.  Flow  from  each  port  was  regulated  with  a rot:u::- 
eter  and  gate  valve  and  could  be  varied  from  0.000^4  to  0.0027  cfs. 
Density  stratification  was  achieved  with  saline  and  fresh  waters.  Con- 
ductivity and  temperature  sensors,  calibrated  against  known  densities, 
were  used  to  measure  density.  Velocity  distributions  were  determined 
from  video  recordings  of  dye  streak  displacement. 

5.  Tiie  test  results  indicated  that,  with  all  four  ports  simul- 
taneously discharging,  the  observed  limits  of  withdrawal  and  release 
densities  could  more  accurately  be  predicted  by  using  the  flow  rate 
and  port  area  for  one  of  the  ports  rather  tiian  tiie  total  flow  rate  and 
port  area.  Tlierefore,  the  Marysville  simulation  model  was  coded  so  tlie 
discharge  through  a single  port,  which  is  obtained  by  dividing  tiie  total 
generation  release  rate  by  6 (for  6 turbines),  is  used  to  compute  the 
limits  of  withdrawal. 
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In  accordance  with  ER  70-2-3,  paragraph  6c(l)(b), 
dated  I5  February  1973>  « facsimile  catalog  card 
In  Library  of  Congress  format  is  reproduced  below. 


